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SUMMARY  REPORT 

This  report  summarizes  the  fallout  pattern  scaling  relationships  that  were 
developed  in  the  period  1962  to  1964;  the  report  includes  the  values  of  the 
scaling  equation  coefficients  that  were  derived  from  selected  fallout  pattern 
data.  The  meaning  of  the  scalar  wind  speed  multiplier  that  is  used  in  the 
scaling  equations  is  discussed  relative  to  computer  applications  of  the  scal¬ 
ing  system  and  approximate  wind  speed  adjustment  factors  for  use  with  wind 
speed  averages  that  may  be  assumed  in  such  applications  are  provided.  The  rel¬ 
ative  degree  of  wind  shear  inherent  in  the  scaling  system  parameters  is  also 
discussed  in  some  detail.  Basic  equations  for  relating  surface  density  of 
radionuclides  and  air  ionization  rates,  Including  consideration  of  fractionation, 
surface  roughness,  and  instrument  response,  are  given  and  discussed  together 
with  the  influence  of  these  factors  and  others  on  the  limiting  values  of  K 
factors  that  represent  the  relative  amount  of  the  radioactive  sources  contained 
within  the  deduced  area  covered  by  the  fallout  patterns.  Scaling  equations  and 
data  are  also  presented  for  use  in  estimating,  for  any  location  in  the  fallout 
region,  the  time  of  fallout  arrival,  the  time  of  fallout  cessation,  the  varia¬ 
tion  of  the  exposure  rate  (i.e.,  all-  ionization  rate  in  roentgens  per  hour) 
with  time  during  fallout  arrival,  and  the  total  exposure  from  the  time  of  fall¬ 
out  arrival  to  selected  later  times. 
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FOREWORD 

The  major  content  of  this  report  was  developed  by  the  author  in  the 
period  1962  to  1964  and,  in  draft  form,  the  material  has  been  available  to 
computer  programmers  at  the  Stanford  Research  Institute  (Menlo  Park,  Califor¬ 
nia)  and  at  the  American  Research  Corporation  (Fullerton,  California).  The 
computer  programs,  in  turn,  have  been  available  to  the  Office  of  Civil  Defense, 
their  contractors,  and  others.  In  some  cases  programming  simplifications  and 
interpolation  schemes  have  been  added  to  decrease  computing  time  or  for  other 
reasons.  Since  such  changes  can  become  iterative  with  respect  to  departures 
from  the  original  systems,  and  since  without  the  original  scaling  functions, 
the  program  user  has  no  means  of  checking  the  program  output,  it  was  requested 
that  the  original  scaling  equations  be  reported  for  record  and  for  computer 
program  verification. 

In  preparing  the  report  for  publication,  a  few  changes  in  scaling  func¬ 
tions  were  made,  mainly  on  the  procedures  for  estimating  the  time  of  arrival 
an-t  rate  of  arrival  of  the  fallout  from  cloud  heights.  Also,  a  few  out-of- 
date  assumptions,  statements,  and  conclusions  were  deleted  or  revised.  New 
work  reported  elsewhere  (such  as  that  sponsored  by  the  Defense  Atomic  Support 
Agency,  Department  of  Defense)  since  1964  is  not  discussed  nor  is  reference 
made  to  such  studies,  since  the  results  therefrom  were  not  available  for  the 
analytical  results  summarized  in  this  report. 

Some  question  still  exists  with  respect  to  the  appropriate  value  of  a 
wind  speed  that  should  be  applied  In  the  calculations  to  conform  with  the  total 
and  angular  shear  effects  included  in  the  derived  scaling  functions.  To  clarify 
this  question,  particle  displacement  calculations  have  been  added,  discussed, 
and  a  first  order  wind  speed  adjustment  factor  is  suggested  for  use  when  an 
average  wind  speed  is  assumed;  the  details  of  the  discussion  on  this  question 
are  given  in  paragraph  3.5.4. 
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Chapter  3 

DISTRIBUTION  OF  LOCAL  FALLOUT 


3 . 1  Background 

3.1,1  The  Distribution  Proceas 


A  very  simple  descriptive  statement  of  the  fallout  process  Is 
that  a  cloud  of  particles  is  formed  rapidly  as  the  result  of  an  explosion 
and  that  this  cloud  is  then  dispersed  by  the  wind  and  by  the  force  of 
gravity  acting  on  the  particles  to  return  them  to  the  earth.  Most  treat¬ 
ments  of  distribution  of  fallout  assume  that  the  visible  volume  occupied 
by  the  nuclear  cloud  and  stem  above  the  point  of  detonation  within  a  few 

minutes  after  explosion  more  or  less  defines  the  volume  source  of  the 
1  2 

fallout  particles.  One  treatment,  however,  considers  the  particle  source 
volume  contained  within  the  air  volume  3wept  through  by  the  rising  fireball. 
In  either  case  the  source  volumes  for  the  particles  depend  on  total  yield 
and,  if  other  than  surface  detonations  are  considered,  on  the  height  or 
depth  of  burst.  The  yield-dependent  parameters  which  are  used  to  define 
the  particle  source  geometry  include  the  cloud  height,  cloud  thickness  and 
radius,  and,  occasionally,  the  stem  geometry,  and  the  time  dependence  of 
these  parameters. 

One  important  additional  factor  that  is  usually  considered  in  the 

fallout  distribution  process  is  the  spatial  concentration  of  the  particles 

in  t.he  volume;  also  qualitative  considerations  have  been  given  to  internal 

3,4 

circulations  of  the  particles  by  several  investigators.  The  discussion 
in  Reference  4  on  this  circulation  is  summarized  in  this  chapter. 

The  fall  trajectory  of  a  particle  through  the  atmosphere  depends  on 

its  own  properties  and  on  meteorological  factors.  The  various  aspects  of 

5  2  1 

these  factors  have  been  discussed  by  Schuert,  Anderson,  and  others. 

The  major  properties  that  influence  a  particle's  fall  rate  through  the 
atmosphere  are  its  density,  diameter  or  size,  and  shape.  The  major  meteo¬ 
rological  factors  are  the  wind  speed  and  direction  and  tne  air  density  and 
viscosity  as  a  function  of  altitude. 
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The  two  air  properties,  of  course,  are  dependent  on  the  air  pressure 
and  temperature,  and  these,  In  turn,  change  with  altitude.  The  wind  speed 
and  direction  arc  also  highly  variable  quantities,  since  each  has  both 
spatial  and  time  variations.  The  vertical  motions  of  the  air  and  particle- 
group  diffusion  can  influence  the  fall  trajectory  of  particles  but  are 
usually  not  taken  into  account  in  the  study  of  the  fallout  distribution 
process . 

It  is  of*en  found  that  when  the  observed  cloud  rates  of  rise  (or  the 
stabilized  cloud  heights)  and  the  particle  fall  rates  are  used  to  compute 
the  time  of  arrival  of  particles  at  locations  very  close  to  ground  jero 
where  fallout  from  stem  altitudes  Bhould  predominate,  the  calculated  arrival 
times  are  quite  consistently  longer  than  the  observed  arrival  times.  Actually, 
the  same  discrepancy  is  often  observed  for  cloud  fallout  at  larger  distances. 

The  consistency  of  the  longer  computed  arrival  time  for  particles  falling 
near  ground  zero  suggests  that  when  the  rising  cloud  takes  on  a  toroidal 
motion,  the  larger  particles  are  thrown  ftom  the  gas  mass  and  experience  down¬ 
ward  accelerations  for  some  rather  extended  period  of  time.  Because  the 
calculated  fall  rates  include  only  accelerations  due  to  gravity,  the  computed 
time  of  fall  (neglecting  downward  accelerations)  from  the  height  of  the  cloud 
would  always  be  longer  than  the  true  time.  Conversely,  when  the  fall  rates 
are  used  in  order  to  estimate  the  height  of  origin  of  a  particle  from  the 
time  of  its  arrival  on  the  ground  (including  its  rise  time),  the  computed 
height  of  origin  i  less  than  the  cloud  height. 

This  interpretation  of  the  above-mentioned  observations  of  particle 
arrival  time  may  be  used  to  describe,  in  qualitative  terms,  the  process  of 
stem  fallout.  The  rising  fireball  takes  on  toroidal  circulation  as  it  rises 
from  the  surface  of  the  ground,  and  this  circulation  persists  through  tran¬ 
sition  to  cloud  form  until  the  internal  pressures  and  temperatures  of  the 
system  approach  those  of  the  ambient  air,  thereby  establishing  a  large-scale 
air  circulation.  Air  and  soil  particles  rise  from  directly  below  the  cloud 
in  a  narrow  visible  stem  or  chimney,  and  the  surrounding  air  is  entrained 
over  the  whole  length  of  this  stem.  This  rising  material  flows  into  the 
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bottom  center  of  the  cloud,  and  the  countercurrent  air  flow,  around  the 
periphery  of  the  cloud,  is  downward.  The  observable  effect,  upon  occasion, 
is  that  the  mass  of  particles  appears  to  flow  out  from  the  top  portion  of 
the  cloud  and  then  downward.  As  the  cloud  approaches  its  maximum  height, 
the  circulation  pattern  apparently  rapidly  disintegrates  or  breaks  up  into 
segregated  regions  of  turbulence  under  influence  of  the  ambient  meteorological 
forces . 

When  the  toroidal  circulation  starts,  a  particle  (or  liquid  drop)  in  the 
central  region  of  the  cloud  would,  by  centrifugal  force,  be  moved  to  the  outer 
periphery  of  the  cloud  and  then  be  accelerated  downward  at  speeds  greater  than 
the  particle's  normal  fall  velocity;  it  would  then  be  at  a  lower  altitude  than 
the  cloud  when  its  terminal  fall  velocity  is  reached.  Even  if  this  centrif¬ 
ugal  action  and  movement  to  the  exterior  of  the  rising  cloud  did  not  occur  for 
the  majority  of  the  particles,  they  could  still  fall  from  lower  altitudes,  by 
virtue  of  the  downward  circulation  around  the  periphery  of  the  cloud,  than 
would  be  calculated  on  the  basis  that  gravity-pull  alone  was  overcoming  the 
gross  rise  r»te  of  the  visible  cloud. 

However,  even  with  toroidal  motion,  the  separation  of  fallout  particles 
by  size  because  of  gravity  forces  is  still  a  valid  concept.  The  smaller 
particles  will  not  move  outward  by  centrifugal  forces  as  far  as  the  larger 
ones  in  the  circulation,  and  they  could  be  swept  back  upward  through  the 
cloud  as  long  as  the  velocity  of  the  rising  all  is  sufficiently  large.  This 
type  of  particle  source  circulation  and  ejection  can  be  used  to  explain  the 
observed  change  In  radiochemical  composition  of  different  size  particles 
discussed  in  Chapter  2,  Reference  6. 

The  major  radiological  factors  in  the  fallout  distribution  process  are 
the  fission  yield  and  the  variation  with  particle  size  of  the  gross  radioac¬ 
tivity  carried  by  particles  of  a  given  size.  The  first  essentially  determines 
the  total  radioactivity  available  for  distribution  on  the  particles;  the  second 
involves  the  distribution  of  that  radioactivity  among  particles  of  different 
sizes , 
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3.1.2  Mathematical  Representations 

The  original  attempt  to  describe  and/or  predict  the  end  result 
of  the  fallou^  distribution  process  —  the  fallout  pattern  —  was  made  by 
C.  F.  K sand a  and  coworkers  in  1953. 7  The  original  scaling  method  was  based 

g 

on  the  work  of  Laurino  and  Poppoff,  which  described  some  fallout  patterns 

for  low-yield  devices  derived  from  observed  data  obtained  during  Operation 

Jangle  in  1951,  The  original  scaling  method  was  intended  for  predictions 

or  estimates  of  fallout  patterns  from  yields  possibly  as  high  as  10  KT. 

9 

In  1955  che  method  was  expanded  to  Include  yields  in  the  megaton  range. 

without  adequate  explicit  experimental  documentation.  This  method  was  sub- 

10 

sequently  included  in  ENW  ;  however,  in  the  latest  edition  of  this  document, 
the  fallout  pattern  scaling  is  revised.  In  many  damage  assessment  studies 
of  fallout  effects,  a  scaling  system  is  to  be  preferred  over  a  complex 
mathematical  model. 

Mathematical  models  attempt  to  establish  quantitative  values  for  the 
several  fallout  distribution  parameters  mentioned  above  and  to  compute  the 
activity  deposited  on  the  ground  at  various  locations,  usually  through  the 
use  of  electronic  computers.  The  general  approach  used  and  the  organizations 
and  investigators  involved  in  the  development  and  testing  of  these  models  up 

to  1957  is  described  In  some  detail  by  Kellogg.1  Later  developments  include 

211  1 2  13 

the  work  by  Anderson,  '  by  Pugh  and  Galiano,  by  Callahan  et  al.,  and 

14 

by  Rapp,  to  mention  a  few  of  the  unclassified  repo* ted  studies.  A  general 
comment  on  the  results  might  be  that  none  of  the  models  agree  with  each 
other  in  several  details  and  that  none  of  the  models  reproduce  very  accurately 
all  of  the  few  data  in  the  yield  range  of  1  KT  to  15  MT  that  are  experimentally 
available . 

The  exact  causes  of  the  differences  among  the  various  models  are  difficult 
to  isolate  for  at  least  two  reasons: 

1.  Each  model  is  different  from  any  other  in  several  of  its 
assumptions  about  para  eter  values  or  in  its  manner  of 
handling  the  many  variables  mathematically. 
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2.  Generally,  the  reports  describing  the  models  do  not 
Include  sufficient  detailed  information  regarding  the 
minor  assumptions  and  the  methods  used  in  making  the 
computations , 

If  the  input  data  in  the  mathematical  models  were  all  more  reliably  estab¬ 
lished  experimentally,  many  of  the  differences  among  them  would  disappear. 
Whether  this  would  produce  better  agreement  with  observations  would  still 
have  to  be  established. 

In  mos  mathematical  models,  after  selection  of  the  values  of  the  source 
geometry,  trajectory,  and  radiological  factors,  the  computation  is  carried 
out  by  dividing  the  source  geometry  for  each  of  several  particle  size  ranges 
into  horizontal  discs  of  finite  thickness.  The  location  on  the  ground  where 
these  "particle  discs"  land,  under  the  influence  of  stated  meteorological 
conditions,  is  then  calculated.  All  the  activity  at  each  of  a  series  of 
coordinate  points  is  then  summed  according  to  the  number  of  different  discs 
that  land  at  the  point  and  the  mount  of  activity  assigned  to  each  disc. 

So  short  a  summary  of  the  work  on  the  mathematical  models  should  not  be 
interpreted  to  mean  that  the  efforts  in  the  development  of  the  mathematical 
models  have  been  small  and  unfruitful.  On  the  contrary,  much  has  been 
learned  about  the  fallout  process  through  them,  and  most  of  the  concepts 
employed  by  many  of  the  mathematical  model  developments  are  covered  in  the 
remainder  of  this  chapter.  But  to  describe  all  the  work  and  all  the  details 
of  each  model  currently  in  use  is  not  considered  to  be  within  the  scope  of 
this  discussion. 

3.1.3  Fallout  Pattern  Features 


Although  observed  data  on  fallout  patterns  f^om  land-surface  shots 
of  various  yields  are  very  meager,  the  processed  data  give  indication  in  a 
qualitatitive  way  of  a  number  of  persistent  characteristics.  For  example,  the 
general  shape  of  the  fallout  standard  intensity  contours  (in  R/hr  at  1  hr) 
from  shots  in  which  the  wind  structures  were  rather  simple,  resembles  a  shadow 
of  the  mushroom  cloud  and  its  stem  on  the  ground. 
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Because  of  the  shortage  of  reliable  data  on  the  fallout  patterns  from 
land-surface  detonations,  any  systematic  method  for  scaling  fallout  patterns 
(l.e.,  methods  for  Interpolating  and/or  extrapolating  data  from  one  weapon 
yield  to  another)  must  take  full  advantage  of  nil  such  apparently  persistent 
qualitative  characteristics  of  the  available  patterns.  In  devising  methods 
that  can  convert  the  qualitative  characteristics  to  quantitative  ones,  the 
methods  must,  of  course,  be  capable  of  at  least  reproducing  the  observed  data 
used  in  obtaining  the  original  scaling  relationships,  which  are  given  as 
functions  of  weapon  yield. 


Some  of  the  apparently  persistent  characteristics  of  the  fallout  patterns 
from  surface  detonations  are: 

1.  In  the  region  near  ground  zero,  the  intensity  gradient  in  the 
upwind  and  crosswind  directions  is  very  steep. 

2.  The  high  intensities  near  ground  zero  appear  as  an  intensity 
ridge  (rather  than  as  a  circular  peak)  displaced  in  the  down¬ 
wind  direction. 

3.  The  length  of  this  high  intensity  ridge  appears  to  be  propor¬ 
tional  to  the  width  of  the  lower  portion  of  the  stem. 

4.  The  peak  intensity  of  the  ridge  increases  with  yield  in  the 
1-  to  10-KT  yield  range  and  decreases  in  the  100-KT  to  10-MT 
yield  range. 

5.  The  best  simple  empirical  relationship  for  the  variation  of 
the  intensity  with  upwind  and  crosswind  distance  from  ground 
zero,  from  graphical  plots  of  the  data,  is  that  the  form  I0e~kx, 
where  IQ  is  the  ridge  peak  intensity,  k  is  a  constant  for 

a  given  yield,  and  x  is  the  upwind  and/or  crosswind  distance 
from  the  upwind  shoulder  of  the  ridge  peak. 

6.  The  contours  downwind  from  ground  zero  appear  to  be  parallel 
to  the  intensity  ridge  for  its  entire  length. 

7.  At  distances  greater  than  the  length  of  the  ridge,  the 
intensity  contours  directly  downwind  decrease  with  distance 
from  ground  zero. 

8.  At  Borne  distance  downwind  (or  perhaps  even  upwind  for  very 
large  yields),  the  low-valued  intensity  contours  fan  out,  and 
the  intensities  directly  downwind  from  ground  zero  rise  sharply 
with  distance  and  then  more  slowly  with  distance  to  a  peak 
value . 
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9.  The  distance  from  ground  zero  to  this  downwind  peak  intensity 
increases  with  weapon  yield. 

10.  The  magnitude  of  the  peak  intensity  also  appears  to  increase 
continuously  with  yield. 

11.  The  distance  between  the  lower  valued  contours  appears  to  be 
related  to  the  width  of  the  cloud  (not  considering  wind  shear 
differences),  and  the  maximum  width  seems  to  occur  farther 
downwind  than  the  peak  intensity. 

12.  The  variation  of  the  intensity  with  downwind  distance  from 
this  outer  pattern  peak  can  be  approximated  within  reasonable 
limits  of  error,  if  the  form  of  the  wind  shear  pattern  is  not 
complex,  by  a  function  of  the  form  I0e-mx,  where  IQ  is  the 
peak  intensity,  m  is  a  yield-dependent  parameter,  and  x  is 
the  downwind  distance  from  the  peak. 


The  above-listed  fallout  pattern  characteristics  are  based  on  a  combina¬ 
tion  of  experimental  observations  and  analyses  of  field  test  data.  The  most 
reliable  persistencies  appear  to  be  those  numbered  1,  2,  4,  6,  7,  8,  and  9. 

3.1.4  Simplified  Fallout  Scaling  System  for  Land-Surface  Detonations 

The  fallout  scaling  system  described  here  was  developed  for 
estimating  standard  intensities,  potential  exposure  doses,  and  other  radio¬ 
logical  quantities  by  use  of  both  manual  and  machine  computational  techniques. 
The  system  is  based  on  corrected  experimental  data,  on  empirical  relationships 
among  the  geometrical  arrangement  of  the  cloud  and  stem  as  the  source  of  fall¬ 
out  particles,  and  on  several  of  the  observed  features  of  the  fallout  pattern 
of  radiation  intensities  on  the  ground.  In  the  system,  the  cloud  and  stem 
dimensions  are  stylized  as  simple  solid  geometric  configurations  to  facilitate 
the  use  of  algebraic  relationships  among  the  model  parameters  and  the  dependence 
of  the  parameter  values  on  weapon  yield. 


In  making  estimates  of  the  hazards  from  fallout,  for  the  purpose  of 
establishing  the  nature  and  required  degree  of  protection  against  these  hazards, 
two  major  quantities  requiring  evaluation  are  (1)  the  exposure  dose  levels 
that  can  result  at  different  distances  from  the  detonation  and  (2)  the  land 
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surface  area  In  which  the  exposure  dose  is  greater  than  a  stated  amount.  To 
make  these  evaluations  requires  estimates  of  the  amount  of  fallout  that 
deposits  at  various  locations,  the  time  at  which  the  fallout  arrives,  and 
the  rate  of  its  arrival. 

Such  general  evaluations  of  radiological  hazard  levels  and  of  the  pro¬ 
tection  requirements  for  radiological  countermeasures  in  defense  planning 
must  first  consider  the  possible  levels  of  effect  (or  hazard)  and,  in  a 
generalized  manner,  the  feasibility  of  methods  for  protecting  against  these 
levels  of  possible  hazard.  For  these  purposes,  a  rather  simplified  fallout 
scaling  system  can  serve;  no  precise  or  accurate  prediction  of  fallout  under 
specified  detonation  and  wind  conditions  is  needed  or  is  possible,  even  with 
the  most  complicated  fallout  models  at  their  present  stage  of  development. 
Therefore  in  the  following  discussion  the  presentation  is  limited  to  the 
description  of  a  simplified  version  of  the  fallout  distribution  process. 

The  mathematical  derivations  of  the  simplified  fallout  scaling  system 
are  designed  to  describe  the  fall  of  particles  of  different  size-groups 
from  a  volume  source  in  the  air;  the  boundaries  of  that  source  are  assumed 
to  depend  only  on  weapon  yield.  The  problem  is  to  describe  mathematically 
the  dependence  of  the  fallout  pattern  features,  in  space  and  time,  on  (a) 
the  cloud  and  stem  geometry,  (b)  the  particle  fall  velocities,  (c)  the  wind 
velocity,  (d)  the  radioactivity-particle  size  distributions,  and  (e)  the 
weapon  yield. 

The  geometrical  configuration  of  the  cloud  for  the  scaling  system  is 
an  oblate  Bpherold,  and  the  configuration  of  the  stem  is  a  frustum  of  an 
exponential  cone  or  horn  whose  larger  base  is  approximately  adjacent  to  the 
bottom  of  the  spheroid.  The  fall  of  particles  from  each  of  these  source 
volumes  is  considered  separately.  The  mathematical  description  and  detailed 
assumptions  used  in  the  development  of  the  model  are  given  in  Reference  4; 
only  those  details  n<  euc  J  for  use  of  the  derived  pattern  scaling  functions 
are  presented  in  the  fallowing  discussion. 
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Some  of  the  pattern  features  of  interest  along  the  downwind  axis  (Y  =  0) 

! 

of  the  Idealized  fallout  pattern  are  shown  In  Figure  3.1  as  a  schematic  inten¬ 
sity  profile.  The  numbers  shown  in  the  figure  correspond  to  the  numerical 
subscripts  of  the  scaling  functions.  The  evaluated  scaling  functions  ;cr 
these  and  other  quantities  are  given  in  the  following  paragraphs  as  summarized 
from  the  data  given  in  Reference  4.  In  the  scaling  system,  the  fallout  pat¬ 
terns  for  the  particles  falling  from  stem  altitudes  and  from  cloud  altitudes 
are  computed  separately;  at  locations  where  these  two  patterns  overlap,  the 
computed  standard  intensities  are  then  added  together.  ThiB  is  illustrated 

in  Figure  3.1  by  the  dotted  line  between  X  and  X_.  The  standard  or  reference 

j  b 

condition  for  all  the  fallout  pattern  scaling  functions  presented  in  the 
following  paragraphs  and  sections  is  100  percent  fission  yield;  the  radioactive 
components  assumed  to  be  present  are  discussed  in  Paragraph  3.4.3. 

3.2  Fallout  Deposition  from  Stem  Altitudes 

3.2.1  Ground  Zero  Intensity  Ridge  j 

The  ground  zero  intensity  ridge  is  depicted  in  Figure  3.1  by  the 

downwind  distances  X  and  X  and  the  intensity  I  ,  The  dependence  of  X  j 

*  J  ^  j  «3  ^  «  ’ 

and  Xg  on  the  average  wind  speed  and  weapon  yield  is  represented  by  ] 

j 

0  230  r  0  1281  5  ■ 

X2  =“  0.0327W  |vw  -  3.96W  Julies ;  W  =  30  to  10  KT  (3.1)  j 

and  -j 

0  230 r  0  5  V 

X3  “  0.0327W  [vw  +  3.96W  Julies;  W  =  30  to  10  KT  (3.2)  j 

for  v  in  miles  per  hour  and  W  in  kllotons  total  yield.  The  values  of  I.  ] 

are  estimated  from  j 


*  Approximate  equalities  Indicate  smoothed  logarithmic  scaling  functions;  equal 
signs  are  used  for  the  basic  scaling  functions  assumed  in  the  derivations. 
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Fig.  3.1.  Schematic  Outline  for  the  Intensity  Profile  Downwind  from  Ground  Zero  Along  the 
Axis  of  the  Fallout  Pattern  (Y  =  0)  as  used  In  the  Simplified  Fallout  Pattern 
Scaling  System 
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_  =  3.02  X  106v  Ww-0-519  R/hr  at  1  hr;  W  =  30  to  105  KT  (3.3) 

2  ,  J  w 

where  I ^  is  the  standard  intensity  for  a  100  percent  fission  yield  weapon, 
and  n^  is  an  empirically  derived  yield-dependent  parameter  whose  value  is 
calculated  from 

=  0.821  -  0.0191  log  W;  W  =  30  to  9,000  KT  (3.4) 


or 


n  -  1.200  -  0.115  log  W;  W  -  9  X  ID3  to  10*  KT  (3,5) 

w 

The  effective  fall  velocity  vector  for  the  median-diameter  particles 
landing  at  (X^  +  X.j)/2  is  given  by 

v  =  ea.3w"0,080  ft/see;  W  =  1  to  105  KT  (3.6) 

2,3 

The  apparent  height  of  origin  of  the  particles  with  the  fall  velocity  vector 
given  by  Equation  3.6  can  be  estimated  from 

z  *  8.04  X  10V-150  feet;  W  =  30  to  105  KT  (3.7) 

i  f  vJ 

3.2.2  Downwind  Intensity  Contours 

The  downwind  distance  to  the  location  X„  of  Figure  3.1  on  the 

4 

center  line  of  the  stem  pattern  is  given  by 

0  203  5 

X.,  =“  0.316v  w  '  miles;  W  =  30  to  10  KT  (3.8) 

4  w 

The  standard  intensity,  for  a  100  percent  fission  yield  weapon,  associated 
with  the  downwind  distance  designated  by  X^  is  given  by 

I.  =  15.0/v  (3.9) 

4  w 
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The  downwind  distances,  ,  to  other  Intensity  values,  1^,  between  1^  ^ 
1^  on  the  pattern  centerline  from  to  X^  are  estimated  from 


and 


X 

c 


<X4  -  X3>  log  (I2  3/Is) 
iog  (I2>3/I4) 


(3.10) 


The  effective  fall  velocity  vector  of  the  media-diameter  particles 

landing  at  X,  Is  given  by 
4 

—0  0ft2  S 

v  =  13. 8W  ‘uo*  ft/sec;  W  =  1  to  10  KT  (3.11) 

4 

The  apparent  height  of  origin  of  the  particles  with  the  fall  velocity 
vector  given  by  Equation  3.11  can  be  estimated  from 

4  0  121  4 

z4  1.57  x  10  W  feet;  W  =  30  to  4  X  10  KT  (3.12) 


or 


z4  ~  1.62  X  loV'118  feet  ;  W  =  4  X  104  to  105  KT  (3.13) 

3.2.3  Upwind  Distance  to  the  1  R/hr  at  1  hr  Contour 

The  upwind  distance  to  the  1  R/hr  at  1  hr  contour  (i.e.,  1^), 
designed  by  the  numeral  one  in  Figure  3.1  is  estimated  from 


0  319  5 

Xx  =  -0.695W  ’  miles;  W  =  30  to  10  KT 


(3.14) 


for  an  average  15-raph  wind  speed.  For  other  average  wind  speeds,  this  distance 
on  the  pattern  center  line  is  calculated  from 

n  e 

X  =  X  -  0.174W  •  log  I  miles;  W  =  30  to  10  KT  (3.15) 

1  a  a  |  J 


where  X  and  I  _  are  given  by  Equation  3.1  and  Equation  3.3,  respectively. 

2,  2  t  o 
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The  distances  to  locations 
other  than  those  of  1^  and  1^  ^ 


between  and  for  the  Intensity  values, 
are  estimated  from 


X' 

c 


X1  + 


tj)  log 


1  /l, 

s  1 


l0«  V^l 


(3.16) 


where  X'  Is  the  distance  from  ground  zero  to  the  Intensity  value  I  along  the 
c  s 

center  line  of  the  pattern. 

Because  the  upwind  fallout  Intensity  contour  locations  are  based  only  on 
empirical  data,  the  scaling  system  cannot  be  utilized  to  estimate  the  diameters 
of  the  particles  landing  at  locations  between  and  X^ . 

3.2.4  Stem  Pattern  Half-Width 

The  (maximum)  half-width  of  the  stem  fallout  pattern,  designated 

by  Ys>  is  the  lateral  distance  from  the  center  line  of  the  stem  pattern  to 

tnc  1  R/hr  at  1  hr  contour;  It  is  located  at  the  downwind  distance  X  given 

s 

by 


X 

s 


V 


V<x2  -  v2  -  v*: 


(X  -  X  )  >  Y 
2  1  s 


(3.17) 


or 


Xs  =  V  <X2  -  V  *  Ys  (3-18) 

The  value  of  Yg(15),  for  an  average  15-mph  wind  speed  is  estimated  from 

Y  (15)  =  0.316W0,400  miles;  W  =  1  to  105  KT  (3.19) 

s 

For  other  wind  speeds,  the  stem  pattern  half-width,  Yg(v^),  is  calculated 
from 
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Y  (v  )  =  Y  (15)S  (v  )  (3.20) 

s  w  s  s  w 

In  which  S  (v  )  is  a  wind  velocity  correction  factor.  The  values  of  S  (v  ) 
s  w  s  w 

are  given  as  a  function  of  weapon  yield  for  several  selected  average  wind 
speeds  in  Figure  3.2,  the  value  of  S^fv^)  for  wind  speeds  other  than  those 
given  in  the  figure  can  be  obtained  graphically  from  a  linear  plot  of  S^v^) 
against  wind  speed  for  a  given  weapon  yield. 

The  lateral  distances  from  the  pattern  center  line  to  other  intensity 

values  between  I  and  the  1  R/hr  at  1  hr  contour  (100  percent  fission  yield) 
2 , 3 

at  the  downwind  distance  X  are  estimated  from 

s 


o 


w  io«  (i2,3/is) 

log  i2  yi1 


(3.21) 


3 .3  Fallout  Deposition  from  Cloud  Altitudes 

3.3.1  Downwind  Intensity  Contours 

The  downwind  intensity  contour  locations  on  the  center  line  of  the 

fallout  pattern  produced  by  particles  falling  from  cloud  heights  are  defined 

in  terms  of  the  intensities  designated  by  I,.,  Ig,  1^,  and  Ig  at  the  downwind 

distances  X_,  X_ ,  X_ ,  and  X  ,  respectively,  as  shown  in  Figure  3.1. 

5  o  7  y 

The  dimensions  and  height  of  the  assumed  elliptical  source  volume  for  the 
particles  falling  from  the  cloud  starting  at  times  of  about  3  to  10  minutes 


after  detonation,  as  derived  from  the  data  reported  by  Pugh  and  Galiano  "  and 
by  Schuert,^  are  represented  by  the  following  yield-dependent  functions: 

r  =  0.464W0-431  miles;  W  =  1  to  105  KT  (3.22) 

b  =  0.265W0'300  miles;  W  =  1  to  105  KT  (3.23) 

h  =  3 . 1 8W° ' 164  miles;  W  =  30  to  105  KT  (3.24) 
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where  r  is  the  cloud  radius,  b  is  the  cloud  half-thickness,  and  h  is  the 
altitude  of  the  center  of  the  cloud.  To  a  first  approximation,  r  represents 
the  cloud  radius  at  the  time  the  cloud  top  reaches  its  maximum  height;  the 
visible  cloud  radius  continues  to  Increase  for  some  time  thereafter. 


In  the  fallout  pattern  scaling  system,  the  variation  of  the  downwind 
distance  with  weapon  yield  of  the  selected  contour  locations  is  mainly 
accounted  for  by  setting  the  various  characteristic  distances  proportional 
to  the  cloud  height  for  a  given  average  wind  speed .  For  the  locations  under 
the  cloud,  some  correction  to  allow  for  a  decreased  effective  cloud  radius 
is  included  in  the  scaling  function.  An  additional,  but  smaller,  dependence 
on  weapon  yield  of  the  distances  to  the  selected  contour  locations  is  included 
to  account  for  a  gradual  shift  in  the  radioactive  content  and  in  the  average 
diameter  of  the  particles  that  land  at  the  designated  locations.  This  change 
in  radioactive  content  is  apparently  caused  by  a  gradual  increase  in  the 
propoition  of  the  radioactive  nuclides  carried  by  the  smaller  particles  as 
the  weapon  yield  Increases. 


The  scaling  functions  for  estimating  the  distances  and  the  intensities 
for  the  selected  contour  locations  are  given  in  two  sets,  depending  on 


whether  the  downwind  distance  is  less 
average  wind  speeds  less  than  a  given 
than  the  cloud  radius,  r,  the  scaling 
speed  limits  for  a  giver,  weapon  yield 


or  greater  than  the  cloud  radius.  For 

speed,  the  distances  X.  and  X_  are  less 

b 

functions  for  the3e  distances  and  wind 
are  given  by 


X  =  0.186W°‘254(v  -  2.49W0,177)  miles;  v  £  3.28W0,245  mi/hr; 

5  w  w 


W  =  30  to  10  KT 


(3.25) 


and 


X  =  0.179W°,309<v  -  2.S9W0’122)  miles;  v  £  2.04W0,231  mi/hr; 

6  w  w 


W  =  30  to  10  KT 


(3.26) 
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The  downwind  distances  to  all  selected  locations  that  are  greater  than 
cloud  radius  are  represented  by 

X  =  v  X®Wni  miles  (3.27) 

1  w  1 

o 

in  which  X^  and  n^  are  constants,  and  the  subscript  i  represents  one  of  the 
selected  locations;  derived  values  of  X°  and  n^  are  summarized  in  Table  3.1. 

The  standard  intensity  for  each  of  the  selected  locations  is  calculated 

from 


I .  =  K°Wrai  log  CD  (v  )  R/hr  at  1  hr 
i  i  i  w 

in  which  K°  and  are  constants,  and  cp^Cv^)  is  given  by 


v  +jv2  +  3.06  v2(m)W0,262 
q  y _ w  7  w  i 

0.0531  v  (m)W°'080[v  +  Jv2  +  1.085W0*102' 

1  w  7  w 

V  s  V  r/h  ;  W  =  30  to  105  KT 
w  i 


(3.28) 


(3.29) 


(v  +  v  r/h) +%v2r2b  2  +  (v  +  v  r/h)2 
,  „  w  i  7  i  w  i 

W  - - r  "  ~  "  ~  ‘  : 

1  w  I  2  2  —2  2 

(v  -  v .  r/h) +%v*r  b  +  (v  -  v  r/h) 
w  i  7  i  w  i 

v  S  v  r/h  ;  W  =  30  to  105  KT  (3.30) 

w  i 

in  wnich  v^  is  the  average  value  of  the  fall  vector  for  the  particles  deposited 
at  the  location  designated  by  1,  and  v^(m)  is  the  minimum  value  of  the  fall 
vector  of  the  deposited  particles.  The  value  of  Ig  obtained  from  Equations 
3.27  and  3.28  results  in 
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Table  3.1 

SUMMARY  OF  DOWNWIND  DISTANCE  AND  STANDARD  INTENSITY  SCALING  EQUATION  CONSTANTS 
FOR  SELECTED  LOCATIONS  ON  TOE  CLOUD  FALLOUT  PATTERN  CENTER  LINEa 


. . . 
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I 


9 


15.0/v 

w 


(3 .31) 


for  v  In  miles  per  hour. 

w 

The  empirically  derived  values  for  the  constants  of  the  intensity  scaling 
functions  for  each  of  the  selected  contour  locations  are  also  summarized  in 
Table  3,1.  The  median  diameter  of  the  particles  deposited  at  each  of  the 
selected  locations  is  estimated  from  the  v  values  calculated  from  the  func¬ 
tions  of  Table  3.1  (after  multiplying  by  1.467  to  convert  the  values  from 
mi/hr  to  ft/sec)  using  data  given  in  Reference  4,  along  with  h  the  cloud 
center  height,  as  the  height  of  origin  of  the  particle  source. 

One  of  the  major  characteristics  of  the  fallout  patter.;  scaling  system 
for  the  fallout  from  the  cloud  altitudes  is  that  it  specifies  a  peak  in  the 
downwind  Intensity  profile  at  X^.  The  intensity,  Iy,  gradually  increases 
with  weapon  yield  approximately  proportional  to  W0*4.  No  experimental  data 
are  available  for  testing  the  reliability  of  the  estimates  of  I?  obtained 
from  extrapolation  of  the  scaling  functions  to  very  high  weapon  yields. 

The  downwind  distances  to  contours  of  other  standard  intensities  on  the 
cloud  fallout  pattern  center  line  are  estimated  from 


(X.  -  X  )  log  (I  /IP) 

X'  =  X  +  -  5  /t  °  I  5  T 

c  5  log  (Ig/I5)  s  6 


(3.32) 


(X,  -  Xc)  log  (I  /Ic> 

V.  _  V  +  7  6 _ 3  6  ■ 

c  *6  +  log  <I7/Ig) 


£  I  s  I_ 
6  s  7 


(3.33) 


cr 


(X  -  X  )  log  (I/I  ) 

X  =  X,  +  —2 - J. - -  ■  -  .  1 — 2-;  i  s  I 

c  7  log  (I?/Ig)  s  7 


(3.34) 


where  X'  is  the  downwind  distance  to  the  intensity,  I  ,  for  the  distances  less 

c  s 

than  X?,  and  Xc  is  the  distance  to  the  intensity,  lg ,  for  distances  beyond  X?. 
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3.3.2  Maximum  Pattern  Half-Width 


The  maximum  pattern  half-width  for  fallout  from  cloud  altitudes 
is  designated  as  Y  ,  and  the  downwind  distance  to  Y  Is  designated  as  X 

o  8  8 

(see  Figure  3.1).  The  empirically  derived  scaling  function  for  X  is  given 

8 

by 


Xa  =  0.325v  w0,315  miles;  W  =  30  to  105  KT  (3.35) 

8  w 

The  crosswind  distances  to  given  contours  in  the  fallout  area  depend, 
first,  on  the  lateral  displacement  of  the  particles  during  the  rise  of  the 
cloud;  second,  on  the  wind  directions  at  all  altitudes  from  the  bottom  to 
the  top  of  the  cloud;  and  third,  on  the  wind  speeds. 

The  wind  speed  has  two  effects  on  the  lateral  displacement  of  an  inten¬ 
sity  contour.  One  is  the  horizontal  displacement  of  particles  with  wind  speed 
because  of  the  relative  horizontal  distance  traveled  in  a  given  period  of  time. 
The  other  is  the  decrease  in  surface  density  of  a  given  size  group  with  wind 
speed  because  of  the  change  in  the  angle  of  the  particle  trajectory.  Hence, 
even  for  the  case  in  which  the  wind  direction  is  the  same  at  all  altitudes, 
a  change  in  wind  speed  results  in  a  change  in  the  maximum  crosswind  distance 
of  a  given  intensity  contour. 

The  values  of  Y°(15)  for  the  maximum  lateral  distance  from  the  pattern 
© 

center  line  to  the  1  R/hr  at  1  hr  contour  for  an  average  wind  speed  of  15 
miles  per  hour,  as  derived  from  both  observed  data  and  ruminations  of  the 
fraction  of  the  radioactivity  contained  in  the  fallout  pattern  (see  Section 
3.4),  are  summarized  in  Table  3.2.  The  computed  pattern  widths  include  the 
effect  of  lateral  wind  shear  contained  in  the  original  data;  this  effect  is 

Q 

discussed  in  Paragraph  3.5.4.  Approximate  scaling  functions  for  Yg(15)  are 
as  follows: 


Y°(15)  =“  0.518W0,615  miles;  W  =  30  to  750  KT 

o 


(3.36) 
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Table  3.2 

SUMMARY  OF  CALCULATED  Y°(15)  VALUES  FOR  SELECTED  WEAPON  YIELDS* 


w 

Y°(15) 

W 

Y°<15) 

(KT) 

(miles) 

(KT) 

(miles) 

5 

1.90 

1 

33.6 

10 

2.28 

2 

40.9 

20 

3.18 

5 

53.4 

50 

5.76 

10 

64.7 

100 

9.10 

20 

78.0 

200 

14.0 

50 

101 

500 

23.5 

100 

123 

a  For  100  percent  fission  yield 
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and 

Y°(15)  =“  4.75W0-283  miles;  W  =  750  to  105  KT  (3.37) 

O 

The  variation  of  Y  with  wind  speed  (for  a  given  wind  direction)  is  deter- 

O 

mined  relative  to  Y„(15)  for  a  wind  speed  of  15  miles  per  hour.  The  represents- 

O 

tion  for  the  variation  is 

W  =  Y8(15)S<Vw)  (3.38) 

in  which  S(v  )  is  the  relative  shear  factor  due  to  wind  speed  only, 
w 

The  values  of  S(v  )  determined  from  the  fallout  scaling  system  parameters 
w 

for  different  wind  speeds  are  essentially  independent  of  weapon  yield .  The 
indicated  value  of  Ya(v  )  is  for  the  particle  groups  falling  at  the  downwind 
distance,  X  ;  the  associated  intensity  contour  that  passes  through  the  loca- 

O 

tion  of  Yg(v^),  Xgt  ls  equal  to  Ig  for  the  same  wind  speed.  The  intensity 

at  the  location  is  thus  1  R/hr  at  1  hr  when  the  wind  speed  is  15  miles  per 

hour.  Values  of  S(v  )  at  several  wind  speeds,  and  the  associated  intensities, 
w 

are  given  in  Table  3.3;  the  tabulated  values  of  s(vw)  are  represented  approx¬ 
imately  by 

S(v  )  *•  0.360(1  +  26.7/v  ) ;  v  =  10  to  22.6  mi/hr  (3.39) 

w  w  w 

and 

S(v  )  “  0.426(1  +  19.0/v  );  v  =22.6  to  50  mi/hr  (3.40) 

w  w  w 

Combining  Equations  3.39  and  3.40  with  Equations  3.36  and  3.37  gives,  for 

Vv„>’ 

Y  (v  )  ="  0.186W0-615  (1  +  26.7/v  )  miles;  v  =  10  to  22.6  mi/hr  (3.41) 

8  w  w  w 

Y  (v  )  ="  0.221W0-615  (1  +  19.0/v  )  miles;  v  =  22.6  to  50  mi/hr  (3.42) 

8  w  w  w 


22 


URS  702-1 


Table  3.3 

SUMMARY  OF  THE  RELATIVE  WIND  SPEED  SHEAR  FACTOR,  S(v  ), 
FOR  THE  FALLOUT  PATTERN  MAXIMUM  HALF-WIDTH  AND  * 
ASSOCIATED  STANDARD  IONIZATION  RATE 
FOR  SEVERAL  WIND  SPEEDS 


Wind  Speed 
(mi/hr) 

S(v  ) 
w 

X9 

(R/hr  at  1 

10 

1.325 

1.5 

15 

1.000 

1.0 

20 

0.840 

0.75 

25 

0.750 

0.60 

30 

0.695 

0.50 

35 

0.658 

0.428 

40 

0.630 

0.375 

43 

0.605 

0.333 
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for  W  values  between  30  and  750  KT ;  and 
0  2  83 

Yc(v  )  ~  1.71W  *  (1  +  26.7/v  )  miles;  v  =  10  to  22.6  mi/hr  (3.43) 

o  w  w  w 

0  2  83 

Y  (v  )  =“  2.02W  *  (1  +  19.0/v  )  miles;  v  =  22.6  to  50  mi/hr  (3.44) 

8  w  w  w 

5 

tor  W  values  between  750  and  10  KT. 

The  value  of  Y°  for  the  1  R/hr  at  i  hr  contour  for  a  given  wind  speed 

O 

(100  percent  fission  yield)  is  estimated  from 


Y°(l) 


W  l0g  *7 

log  a7/i9) 


(3.45) 


where  I  and  I  are  the  intensities  at  X  and  X  (Y  =  0)  for  a  given  value  of 
7  9  7  9 

v  . 
w 

The  maximum  lateral  distance  to  other  intensity  values  between  I?  and 
the  1  R/hr  at  1  hr  contour  at  X  ,  Y°(l)  is  estimated  from 

o  8 


Y8U)  =  Yg(  1 ) 


log  (I./l  ) 
/  s 

log  Iy 


(3.46) 


The  downwind  distance  to  Y  (I)  is  given  by 

O 


Vx)  =  x7  + 


(X  -  X  )  log  (I_/I  ) 

o  T  f  S 

log  (I7/Ig) 


(3.47) 


3.4  Idealized  Fallout  Pattern  Shapes  and  Areas  Within  Contours 
3.4.1  Stem  Pattern  Contours 


With  the  simplified  fallout  scaling  system  and  the  stylized  down¬ 
wind  intensity  profile,  it  is  convenient  to  construct  contours  by  using  simple 
geometric  forms  which  approximate  to  some  degree  the  true  shapes  of  the  contours 
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and,  within  reasonable  limits,  account  .for  the  fraction  of  the  activity 
produced  that  falls  back  to  earth  within  their  areas. 

The  assumed  shape  of  the  intensity  contours  for  the  fallout  from  stem 

altitudes  is  a  combination  of  a  half-ellipse  and  a  partial  circle  having  its 

center  at  X.,  providing  the  distance  (X_  -  X')  Is  equal  to  or  larger  than  the 
2  2  c 

distance  Yq  for  the  same  contour  of  the  stem  pattern.  For  these  conditions, 

the  major  axis  of  the  half-ellipse  on  the  center  line  in  the  downwind  direc- 

o  a  2  1/2 

tion  is  (X  -  x  )  where  x  is  equal  to  X.  +  [(X.  -  Y')“  -  Y  )  ]  and  the 

C  S  S  <2  6  C  O 

minor  axis  is  Y°  the  ellipse  is  centered  at  on  the  pattern  center  line. 

The  equations  for  the  assumed  intensity  contour  shapes,  accordingly, 
are  given  by 


(X 


V 


Y2  =  <X„ 


-  X’>‘ 

C 


X  -  -X’ 
c 


to  X 


(3.48) 


and 


<*  -  v 

(xc  -  *s>s 


<Y°)2 

o 


=  l;  X 


x  to  X 
s  c 


(3.49) 


The  area  contained  within  a  given  standard  intensity  contour,  for  the  above 
assumed  contour  geometry,  is  given  by 


a  =  3 . 1416(X  -  X')"2  1 

s  2  c 


sin 


■‘fc 


/(X„ 


180 


1 . 5708Y 


-  0.5708Y 


Y2;  (X  -  X’)  2  Y° 
o  2  c  o 


(3.50) 
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or 


l0«  <I2,3/Is) 

3.1416(X2  -  XL)2 

sin_1Y  (v  )/(X„  -  X.) 

,  s  w  2  1 

“S  "  l°*  X2.3 

l0g  *2,3 

180 

1.5708Y  (v  )(X„  -  X0) 
s  w  4  3 

iog  a2yv 


0.5708Y  (v  ) 

9 


log  I 


IV  ;  r 

2-2-  J, 

2,3  ' 


(X  -  XI2  -  Y2(v  ) 
2  1  s  w 


log  (I2,3/Is) 


4-  1.5708Y  (v  )(X,  -  X.) 

s  w  3  2 


(X, 


-  X, )  2  Y  (v  ) 
1  s  w 


(3.51) 


The  total  activity  contained  within  the  stem  fallout  pattern,  obtained 
from  the  integration  of  I^da^  (where  da^  is  calculated  from  Equation  3.51  over 
the  intensity  limits  from  I  to  I  ,  is  given  by 

C  I  o 


A 

s 


*2,3 

1.186(X2  -  Xx)2 

sin""1Y  (v  )/(X0  -  X) 
s  w  2  1 

108  *2,3 

108  *2,3 

1  180 

h~  — 

0.593Y  (v  )(X  -  X_) 

_ s  w  4  3 

iog  <I23/I4) 


0.2154Y  (v  ) 
s  w 


V 


{X2  ' 


V 


-  Y2(v  ) 
s  w 


log  I 


2,3 


+  0.682Ys(vw)(X3  -  X2) 


<X2  -  V  *  W!  *2,3  ^*4 


(3.03) 


For  some  conditions  of  wind  speed  and  weapon  yield,  the  maximum  half¬ 
width  distance,  Y°,  will  exceed  the  distance  (X2  -  X^).  As  an  example, 

Y  (v  )  for  the  1  R/hr  at  1  hr  contour  exceeds  the  distance  (X^  -  X  )  for 
s  w  2  1 

weapon  yields  larger  than  22  WT  and  a  15  mile  per  hour  wind  speed.  The 
assumed  contour  shape  for  the  case  where  Y°  2  (X„  -  X')  is  that  of  two  half- 
ellipses,  both  centered  at  the  point  X,,  on  the  pattern  center  line.  The 
upwind  half-ellipse  has  a  major  axis  given  by  Y°  and  a  minor  axis  given  by 
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(X  -  X');  the  downwind  half-ellipse  has  a  major  axis  given  oy  (X  -  X  ) 

“  o  c  2 

and  a  minor  axis  given  by  Y  .  The  general  equations  for  the  ideali/.cd 
intensity  contour  shapes  are 


<X2  ~  X> 

(x2  -  x;)s 


(Y°)2 

o 


=  l ;  x  £  x„ 


and 


(X  -  X2)2  y2 

(X  -  X  )2  <Y°)2 

c  2  o 


X  *  X„ 


(3.53) 


(3.54) 


The  area  contained  within  a  given  standard  intensity  contour,  for  the 
two  half-ellipses,  is 


a  =  1 . 5708Y (X  -  X*  )  :  Y  *  (X  -  X*  ) 
s  c  c  o  2  c 


(3.55) 


1.5708Ys(vw)  log  d2  3/Is) 

_(X2  -  V 

«4  -  X3>  ' 

%  log  I2  3 

_l0g  l2 ,3 

l0g  ^I2,3/I4^ 

iog  (i2yV 


+  «3  -  V  ;  W  2  (x2  -  V 


(3.56) 


The  total  activity  contained  within  these  contour  areas  is  given  by 


0.682Y  (v  )I 
.  s  w  2,3 

S  -  log  I23 


(X,  -  X  >  (X  -  X  ) 
0.858  j - —  +  4  J 


log  I23  log  (X23/I4)| 


+  (X3  -  X2); 


W  ^  (X2  *  Xl>:  *2,3  ^  *4 


(3.57) 
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3.4.2  Cloud  Pattern  Contours 


The  assumed  shape  of  the  intensity  contours  for  the  fallout  pattern 
formed  by  particles  filling  from  cloud  altitudes  is  described  by  two  half- 
ellipses  smoothly  joined  at.  the  downwind  distance  of  the  maximum  width  of  the 
contour.  The  locations  of  the  maximum  contour  widths  fall  on  a  line  joining 
X  on  the  pattern  center  line  (Y  =  0}  and  the  point  X.(I) ,YQ(I )  (see  Equations 
3.46  and  3,47)  . 

The  general  equations  for  this  idealized  contour  shape  are 


[X  (I)  -  X]2  2 

- -  +  - -  =  1;  X  re  X  (I)  (3.58) 

[X  (I)  -  X^ ]  \2(I) 


and 


[X  -  X  (I)]2  2 

- -  +  -I -  =  1;  X  2  X  <I)  (3.59) 

[X  -  x  u)]2  y2(D  8 


The  area  contained  with  a  given  standard  intensity  contour  is  given  by 


a  =  1.5708Y  (X  -  X')  (3.60) 

c  8  c  c 

in  terms  of  the  designated  distances;  in  terms  of  intensities  and  scaled 
pattern  features,  the  areas  within  contours  are  given  by 

(Xg  -  X7)log(I7/Is)  (X?  -  X6)log(Is/I6) 
log  <I7/Ig>  log 


1.5708Y8(vw)log(I7/Is) 
log  U7/Ig) 


+ 


;  I 

s 


£  I, 


(3.61) 
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and 


a 

c 


1 -5708Y8(v^)log(I7/Is) 
'log  a7/i9) 


-  X7)log(I7/Ig) 
log  a7/i9) 


(x.  -  x  "oga  /i.) 

7  8  P 

lo6  (I7/I9) 


+ 


(3.62) 


The  activity  integral  for 
be  estimated  from 


the  cloud 


fallout  pattern,  designated  A  , 

c 


can 


A  = 
c 


1.364Yfl(v  ) 

8  w 

(X9-V 

log  I? 

log  (I?/Ig) 

7  6 


0.434(I7  -  I6)  -  I6  log  (17/I6) 


°.6822Y8(vw)I6 

(x3  -  X?) 

log  I? 

log  (VV 

0 . 86 8 J 


(X6  *  X-} 


l0g  (V-5 


f-y  [log  (I7/I6)  -  log  <Ig/I5)  +  0 . 868 J  +  X?  -  X5|  (3.63) 


for  the  case  where  »  1^ 


3.4.3  Summary  of  Pattern  Features  and  Activity  Summations 

The  fission-product  radioactivity  produced  by  a  detonation  is  about 
23 

1.4  X  10  BW  fissions  where  W  is  the  total  yield  in  KT  and  B  is  the  ratio  of 
fission  to  total  yield.  Some  fraction  of  this  total  activity  is  contained 
within  ‘he  fallout  pattern.  If  the  activity  in  a  pattern  is  summed  over  the 
fallout  area  from  the  central  high  Intensities  down  to  a  stated  low-valued 
Intensity  contour  enclosing  the  largest  area,  the  fraction  accounted  for 
increases  with  yield.  An  ionization  rate  or  intensity  sum  of  a  fallout  pattern 
made  in  this  way  does  not  account  for  the  activity  deposited  on  the  ground  at 
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lower  Intensities  than  the  selected  "lowest”  contour,  nor  does  it  account  for 
the  fraction  carried  away  on  very  small  particles  as  world-wide  fallout.  The 
sum  or  integration  of  the  activity  over  the  fallout  area  can  be  used,  however, 
to  determine  the  fraction  of  the  activity  that  is  accounted  for  in  the  fallout 
pattern  out  to  a  stated  low-level  contour;  the  equations  of  such  a  summation 
for  the  simplified  fallout  scaling  system  are  given  in  the  two  preceding 
sections . 

The  amount  of  activity  accounted  for  in  a  fallout  pattern  is  defined  by 

a 

_  23  C  ° 

C(1)K(1)  1.4  X  10  BW  =  \  I  da  (3.64) 

Jo  s 


in  which  C(l)  is  the  true  fraction  of  the  number  of  fissions  accounted  for, 
and  K(l)  is  the  average  value  of  the  ratio  of  R/hr  at  1  hr  to  the  number  of 
fissions  or  kiloton  equivalents  per  unit  area  for  all  the  activity  within 
the  area,  aQ.  The  average  or  accumulated  fractionation  number,  for  the 
radioactive  elements  accounted  for,  may  be  defined  through  K(l)  by 


K(  1 )  =  5q  [rfp(l)  ifp(l)  +  ijUjj  (3.65) 

in  which  5  is  an  instrument  response  factor,  q  is  a  terrain  attenuation  factor, 


r^  (1)  is  an  average  gross  fractionation  numbe1*  for  the  fission  products  appli¬ 
cable  to  the  ionization  rates  at  H  +  1,  1^(1)  is  the  H  +  1  air  ionization 
rate  at  3  ft  above  an  infinite  ideal  plane  for  a  unit  surface  density  of 
unfractionated  fission  products  spread  over  the  plane,  and  i^(l)  is  in  the 
same  units  for  neutron-induced  activities.  The  units  of  measure  for  K(l) 

are  the  same  as  for  i.  (1)  and  i,(l)  and  usually  are  given  in  terms  of  either 

fp  i 

R/hr  at  1  hr  per  fission/sq  ft  or  R/hr  at  1  hr  per  KT/sq  mi.  It  may  be  noted 
that  the  true  air  ionization  rate  is  obtained  when  the  value  of  6  is  set  equal 
to  one . 
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The  empirical  constants  for  the  scaling  functions  of  the  simplified  fall¬ 
out  scaling  system  apply  to  the  reference  condition  of  100  percent  fission 
yield  (B  =  1) ,  Including  0.8  neutron  captures  per  fission  by  U-238  to  give 
the  appropriate  value  of  Ig  in  R/hr  at  1  hr. 

Also,  the  I  values  (or  standard  Intensities)  correspond  to  radiac  measure- 
s 

ments  taken  at  3  ft  above  an  extended  open  area  contaminated  with  fallout  as 
it  existed  when  the  measurements  were  taken.  The  reference  radiac  instrument 
for  the  Ig  values  is  the  AN/PDR-39(T1B)  portable  radiac;  it  has  a  geometric 
and  photon  energy  response  at  very  nearly  75  percent  of  the  true  air  ionization 
rate  at  3  ft  above  a  plane  source  of  fission  products  uniformly  distributed  on 
the  area.  Therefore  to  obtain  the  true  air  ionization,  in  calculating  exposure 
doses,  the  Ig  values  should  be  multiplied  by  1.33. 

The  average  value  of  the  terrain  shielding  factor,  which  is  automatically 

contained  in  the  source  data,  is  about  0.75.  The  data,  in  general,  apply  to 

-13 

U-238  fission,  for  which  the  i^  (1)  of  Equation  3.65  is  equal  to  6.94  X  10 
R/hr  at  1  hr  per  fission/sq  ft  or  3,510  R/hr  at  1  hr  per  KT/sq  mi.  The  value 
of  ii(l)  for  the  indicated  induced  activities  is  0.13  x  10  13  R/hr  at  1  hr 
per  fission/sq  ft.  With  these  numbers,  Equation  3.65  becomes 


K(l)  =  3.90  x  10 


13  [ 


r(l)  + 


0.019] 


R/hr  at  1  hr 
fission/sq  ft 


(3.66) 


or 


K<1)  ■  1,970  [r<l)  ♦  0.019]  (3.67) 

With  Equation  3.67,  the  intensity-area  integral  becomes 

Igda  =  1,970  £r(l)  +  0.019^  C(1)BW  (3.68) 

0 

It  can  be  seen  that,  for  unfractionated  fission  products  (r(l)  =  1)  and 
for  C(l)  equal  to  unity,  the  ideal  or  upper  limit  value  of  the  intensity-area 
integral  per  unit  yield  would  be  about  2,010  R/hr  at  1  hr  for  observed  fallout 
pattern  intensities. 
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Values  of  the  various  fallout  pattern  features  from  the  various  scaling 
functions  and  intensity-area  integral  parameters  are  summarized  in  Table  3.4 
for  selected  weapon  yields.  The  area  integrals  indicate  that  the  fraction  of 
the  total  activity  contained  within  the  1  R/hr  at  1  hr  contour  for  the  simpli¬ 
fied  fallout  pattern  scaling  systems  varies  from  about  0.4  for  a  weapon  yield 
of  1  KT  to  about  0.8  for  a  weapon  yield  of  100  WT .  Although  the  scaling 
functions  given  in  the  previous  sections  apply  only  to  weapon  yields  in  the 
range  of  30  KT  to  100  ifT  (for  the  latter  high  yield  as  an  extrapolation), 
tabulated  values  for  1  and  10  KT  are  included  in  Table  3.4  as  derived  from 
the  data  of  Reference  4. 

3.5  Ionization  Rates  and  Potential  Fxposure  Doses  During  Fallout  Arrival 
3.5.1  Times  of  Fallout  Arrival  and  Cessation 


The  arrival  time  of  fallout  particles  from  stem  altitudes  is 
estimated  as  the  sum  of  two  time  periods.  The  first  of  these  is  the  time 
required  for  the  particles  of  a  given  average  fall  velocity  to  pass  through 
the  fireball  and  circulate  in  a  toroidal  motion  until  the  updraft  from  the 
buoyant  gases  no  longer  carry  the  particles  upward;  the  second  time  period 
is  that  of  gravity  fail  under  influence  of  winds  to  the  ground. 

The  apparent  height  at  which  particles  of  different  diameters  (and  given 
fall  velocities)  are  effectively  free  of  the  updraft  under  the  rising  fireball 
and  cloud  is  estimated  from 

Zf  =  1.80  X  loV’116  £l  -  exp  (-0.00599W-0,116t)j  ft; 

t  2  20  sec  ;  W  =>  30  to  105  KT  (3.69) 

The  rate  of  increase  of  the  apparent  height  with  time  after  detonation  at 
which  the  particles  assume  free  fall  is  then  given  by 
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Table  3.4 

SUMMARY  OF  FALLOUT  PATTERN  FEATURES  AND  FALLOUT  SCALING  SYSTEM 
PARAMETER  VALUES  FOR  AN  ASSUMED  EFFECTIVE  WIND  SPEED  OF  15  MFH 


Pattern  Feature  Weapon  Yield 


or  Quantity 

1  KT 

10  KT 

100  KT 

1  MT 

10  MT 

100  MT 

X1 

-0.385 

-1.21 

-3.01 

-6.29 

-13.1 

-27.3 

0.0648 

0.342 

0.735 

0.882 

0.574 

-1.02 

S 

0.352 

0.967 

2.06 

4.00 

7.68 

14.7 

1.46 

5.18 

12.1 

19.3 

30.8 

48,6 

*. 

0.834 

2.58 

5.00 

7.03 

4.39 

-14.3 

1.34 

4.05 

8.56 

13.6 

21.6 

27,9 

*7 

1.38 

5.30 

14.4 

28.8 

58.4 

118 

*3 

1.91 

7.56 

20.8 

42.8 

88.6 

183 

*9 

29.4 

61.2 

126 

258 

528 

1,080 

Ys(15) 

0.316 

0.796 

1.99 

5.02 

12.6 

31.6 

Yg( 15) 

1.26 

2.28 

9.10 

33.6 

64.7 

123 

X2,3 

1,730 

22,900 

35,300 

9,800 

3,550 

1,440 

I6 

130 

222 

540 

1,720 

5,270 

15,500 

I7 

180 

317 

803 

2,712 

9,070 

30,510 

7  b 

A  /W 
s 

93.4 

724 

544 

77.5 

14.4 

3.3 

A  /Wb 

c 

757 

416 

687 

1,250 

1,412 

1,552 

AT/Wb 

850 

1,140 

1,231 

1,328 

1,426 

1,555 

(  r  ( 1 )"  +  0 ,019)C 

0.785 

0.787 

0.789 

0.793 

0.797 

0.803 

C(l) 

0.550 

0.735 

0.792 

0.850 

0.908 

0.983 

d 

Z2,3 

3,430 

8,910 

15,800 

22,900 

32,300 

44,400 

z  6 

5,200 

15,600 

30 , 800 

40,200 

52,500 

68,500 

hd 

6,600 

18,400 

35,700 

52,300 

76,200 

111,000 

a  Distances  are  in  miles;  intensities  are  in  R/hr  at  1  hr  (observed); 
reference:  100  percent  fission  yield 

b  In  R/hr  at  1  hr  per  KT/sq  mi;  Integrated  to  the  1  R/hr  at  1  hr  contour 

c  r(l)  calculated  from  data  of  Chapter  4 

d  Altitudes  in  feet 
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Zf  =  108  exp  (-0.00599w"°,U6t)  ft/sec;  t  *  20  sec; 

W  =  30  to  105  KT  (3.70) 

Approximate  representations  of  the  relationships  between  the  particle  fall 
rate,  Z^,  at  the  height  Z  and  the  particle  fall  velocity  vector,  v^,  for  a 
particle  falling  from  the  height  Z  to  the  ground  (sea  level)  were  derived  for 
spherical  particles  from  particle  fall  rate  data  as  computed  by  Cassidy15; 
these  are; 

*  •  c  c 

Zf/vf  =  1.00  +  2.00 x  10  Z{;  v  <  2.0  ft/sec;  Z  =  0  to  10  ft  (3.71) 

Zf /vf  =  0.948  -  0.0530  log  v{  +  3.78  X  10”6(1.259  +  log  vf)Zf; 

vf  =  2.0  to  200  ft/sec;  Z  =  104  to  6  x  104  ft  (3.72) 

The  relatlonsnlp  among  Zf,  vf,  and  Zf  given  by  Equation  3.72  represents 
the  tabulated  data  within  10  percent  at  the  extremes  of  the  indicated  limits 
of  application  for  vf  and  Z.  The  indicated  ranges  in  these  parameters  include 
all  the  values  applicable  to  the  fallout  from  the  stem  heights  of  interest; 
therefore,  Equation  3.72  generally  is  used,  along  with  Equations  3.69  and 
3.70,  to  estimate  the  time  after  detonation  when  particles  with  a  given  fall 
vector  leave  the  rising  cloud  and  take  on  free  fall.  This  time  is  given  by 

tx  =  0. 107W° • 116  |log[l585/vf  +  (1.259  +  log  vf)  W0,116] 

-  log  [o. 779(17. 887  -  log  vf)  +  (1.259  +  log  vf)  W°-116]|hr  (3.73) 

The  time  for  the  particles  with  a  given  value  of  v  to  fall  from  the  height 
Z ^  (Z^  £  60,000  ft)  is  estimated  from 
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t„  = 


3.895  ^2 , 034/ v^  +  log  vf  -  17.887J 


F 


585W 


-0.116 


+  vf  (1 .259  +  log  vf)] 


hr 


(3.74) 


Equations  3.73  and  3.74  apply  to  v  values  between  2  and  200  ft/sec  and  to  W 

5  * 

values  between  30  and  10  KT.  With  Equation  3.71,  the  times  are 


t,  =  0.1.07W 


0.116 


log  (27.78  + 


W0-116)  - 


log  (2 , 995/ v ^  + 


W 


0.U6,] 


)  hr 


(3.75) 


and 


f-  —  . , 

277.8 

108  -  v 

t _ La 

t2 

Vf 

2.0  vf  + 

5,990.-°-n6] 

(3.76) 


for  v^  values  less  than  2  ft/sec  and  W  values  between  30  and  10J  KT . 

The  time  of  arrival  on  the  ground,  t  ,  for  the  particles  with  the  falling 

a 

velocity  vf  is  equal  to  t^  +  . 

The  particles  that  reach  the  earth's  surface  first  are  those  with  the 

falling  velocity,  v,-  .  For  these  particles  the  times  after  detonation 

12,3 

of  leaving  the  cloud  and  of  fall  from  Zf  to  the  ground  are  given  by 


t~  =  0.107W 


290.1W°‘080  +  (38.66  -  log  W)W0,116] 


°-U6  (log[ 

-  log  [o.77S  (200.7  +  log  IV)  +  (38.66  -  log  W)W°'116]|hr ; 


W  =  30  to  10  KT 


(3.77) 


and 


t°  =  0.0570W0,080 


.372 .2W 


0.080 


-  200.7  -  log  W J 


-0.036 

290. 1W  +  38.66  -  log 


hr;  W  =  30  to  10  KT 


(3.78) 


Values  of  the  earliest  fallout  arrival  on  the  ground  from  stem  heights  as 
calculated  from  Equations  3.77  and  3.78,  are  given  in  Table  3.5  for  several 
weapon  yields. 
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The  calculation  of  the  time  of  fallout  arrival  and  of  fallout  cessation 
(and  the  time  period  of  deposition)  for  a  given  location  on  the  stem  fallout 
pattern  canter  line  requires  estimates  of  the  minimum  and  maximum  values  of 
Vj  that  arrive  at  the  location.  These  two  values  of  v  are  for  the  particles 
that  fall  from  the  upwind  and  downwind  sides  of  the  stem  volume;  they  are 
determined  from  the  stem  radius  and  various  relationships  between  the  down¬ 
wind  distance,  X,  the  fall  vector,  v^,  and  the  apparent  height,  at  which 

the  particles  take  on  free  fall. 


Table  3.5 


ESTIMATED  TIMES  OF  EARLIEST  FALLOUT  ARRIVAL  FROM  STEM  ALTITUDES 


w 

O 

t  , 

t° 

t0 

t° 

1 

2 

a 

a 

(Kr) 

(hr) 

(hr) 

(hr) 

(min) 

102 

0.062 

0.098 

0.160 

10 

103 

0.093 

0.167 

0.260 

16 

104 

0.135 

0.280 

0.415 

25 

10s 

0.192 

0.463 

0.655 

39 

The  maximum  and  minimum  downwind  distances  of  deposition  along  the  center 
line  of  the  stem  fallout  pattern  for  the  particles  that  fall  freely  from  the 
height,  Zf,  are  estimated  from 

4  0  333  /  -4  -0  131  \ 

X  =  2.73  X  10  (v  /v  )Z  ±  0.0395W  *  exp  (l.48  x  10  W  '  Zj  miles; 

w  f  f  \  f  / 

W  -  30  to  105  KT  (3.79) 

where  v  is  in  mi/hr;  v.  is  in  ft/sec;  and  Z,  is  in  feet.  The  value  of  Z,  for 
w  f  f  f 

a  selected  value  of  v  is  estimated  from 
106(108  -  v  ) 

=  - - -  -  ft;  v  <  2  ft/sec;  V  =  30  to  10  KT  (3.80) 

1  (2.0v,  +  5.990W-0-116)  f 
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or 


1.40  X  104^2,034  -  (17.887  -  log  vf)vf] 


¥■ 


259  +  log  vj)v£  +  1.585W 


-0.116 


] 


ft  ; 


v  =  2  to  200  ft/sec;  W  =  30  to  105  KT  (3.81) 

The  minimum  and  maximum  values  of  v  are  determined  by  selecting  a  series 
of  values  of  starting  with  v^(2,3)  as  the  largest  value  and  then  computing 
Z^  and  the  two  values  of  X  in  Equation  3.79;  the  v  values  are  plotted  as  a 
function  of  X  to  determine  the  paired  values  of  v  at  a  given  value  of  X. 

Such  a  plot  is  shown  in  Figure  3.3  for  the  fallout  from  stem  heights  of  a  5-WT 
yield  detonation.  The  time  of  arrival  of  the  particles  with  v^  values  from 
the  minimum  to  the  maximum,  for  a  given  value  of  X,  is  computed  by  substituting 
selected  values  of  v  in  Equations  3.73  and  3.74  or  Equations  3.75  and  3,76 
and  solving  for  both  t^  and  . 

Estimates  of  the  v  values  and  arrival  times  of  the  particles  deposited 
at  locations  which  are  a  distance  of  Y  miles  from  the  center  line  of  the  stem 
fallout  pattern  can  be  made  by  replacing  the  second  term  of  Equation  3.79 
with 

±  J^1.56  X  10-V-667  exp  (2.96  X  10~ 4W_C ' 131Zf )  -  Y2]  (3.82) 

The  particles  falling  from  cloud  altitudes,  the  time  of  arrival  on  the 
ground  of  particles  with  the  falling  velocity,  v  ,  may  be  generally  represented 
by 


t 

a 


0 . 17W~°  * 1  +  t 


hr 


(3.83) 


-O’ 

where  0.17W  hr  is 
time  of  fall  from  the 
below  in  relation  to: 


the  apparent  time  to  cloud  stabilization  and  t  is  the 
height  Z.  The  application  of  Equation  3.83  is  discussed 
(1)  the  estimation  of  the  time  of  arrival  and  cessation 
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Fig.  3.3.  Variation  of  the  Computed  Minimum  and  Maximum  Values  of  vf  with  Downwind  Distance 
(center  line  of  stem  pattern)  for  Particles  Falling  from  Stem  Altitudes 


UKS  702-L 


of  fallout  at  a  given  coordinate  point,  X,  with  an  assumed  reference  pattern 
spread  due  to  lateral  wind  shear;  (2)  the  estimation  of  the  range  in  falling 
speeds  (and  diameters)  of  particles  that  arrive  at  a  given  location  (as  needed 
for  estimating  the  gross  activity-particle  size  distribution);  and  (3)  the 
estimation  of  the  influence  of  lateral  wind  shear  ut  the  cloud  heights  on  the 
spread  of  the  fallout  pattern.  In  all  treatments,  it  is  assumed  that  the  cloud 
radius,  r,  remains  constant  in  the  direction  of  the  wind  (i.e.,  the  continuing 
lateral  expansion  of  the  cloud  after  0.17W  hr  is  neglected). 

The  earliest  arrival  time  of  particles  falling  from  the  so-called  stabi¬ 
lized  cloud  is  that  of  the  largest  particles,  with  a  fall  vector  designated 
by  v,.,  from  the  lower  surface  of  the  cloud.  The  general  area  on  which  these 
particles  fall  is  defined  by  an  ellipse  whose  axis  on  the  center  line  of  the 
pattern  is  equal  to  r  and  whose  axis  in  the  cross-wind  direction  is  equal  to 

a  distance  designated  as  V  .  The  first  arrival  time  at  downwind  locations 

o 

outside  this  elliptically  shaped  area  of  particles  with  falling  speeds  less 
than  v^  is  for  the  particles  falling  from  the  forward  edge  of  the  cloud  'at 
the  cloud  center  height,  h).  The  general  equation  for  of  the  particles 
falling  from  the  lower  surface  of  the  cloud  that  land  at  the  point  X,  V,  is 


t  „  = 


(vfhrV2  +  Vwx) 

i- r v 

1/2  2, 
vV  1 

-2  2\ 
)  +  V 

vj 

|(l  -V 

2  2) 
/Yo) 

|-  (v  h  1x  -  v  'j 
\  i  w/ 

2  2-2 
|  ll  b 

1 

f  2  2-2 
v  r  b 
f 

1 

hr  (3.84) 


in  which 


2  -2  0.426  5 

hr  b  =  9.75W  mi;  W  =  30  to  10  KT 

0.431  5 

r  =  0.464W  mi;  W  -  1  to  10  KT 

2  -2  -0.262  5 

r  b  =  3.07W  mi;  W  =  1  to  10  KT 

h_1  =  0.314W-0,164  mi"1;  W  =  30  to  105  KT 

2  -2  -0.272  5 

h  b  =  144W  ;  W  =  30  to  10  KT 


for  v  and  v,  in  mi/hr. 
w  f 
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The  assumed  equation  for  Y  ,  the  peripheral  limit  of  the  fallout  pattern 
in  the  cross-wind  direction  is 


Y 

o 


2Y„(  v  )  1 

8  w  L 


-k  (X  + 
e  c  c 


mj  les 


(3.85) 


where  Yg(v^)  is  given  by  Equations  3,41  through  3.44;  this  form  of  represen¬ 
tation  assumes  persistence  of  lateral  shear  such  that  the  pattern  spread  at 

very  large  downwind  distances  is  twice  the  wiuth,  Y  (v  ),  at  the  downwind 

8  w 

distance,  X0.  Further,  it  is  assumed  that  Y  is  equal  to  1.1  Y„(v  )  at  X„ 

8  o  8  w  8 

so  that 

k  =  0.8/(X„  +  r)  miles”1  (3.86) 

c  8 


or 


k 

c 


2 . 46W 


-0.315 


(v  +  1.395W 


,0 . 1 16\ 


mi  les 


-1 


) 


(3.87) 


Substitution  of  v^  from  Table  3.1 


for  v 


in  Equation  3.84  gives 


^ISW0*404  ♦  V  X^ 

t  „  \ _ 'LJ- 

f  0,218  2 

1 , 550W  +  v 

w 


0.464W°’43)  l,550W0-°l3-t-  v3^l  -  Y2/Y^  -  144^7,08  X  W  °,186-v^2V» 


-0.272 


0.218  ? 
1 , 550W  +  v 

w 


0  o  a  -  ^ 

hrs  ;  v  3 .  82VV  ;  W  =  30  to  10  KT 

v. 


(3.88) 
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(l67„°'336  *  vwx) 

,f  '  89S.0-082  *  v3 

w 


0.464.0-431  V^6W°-082  4»;)(l  -  V3/,;)-  144(s.38  X.'0-254-  .J1. 


0,272 


0.082  2 
896W  •  +  v 

w 


0  245  5 

hrs ;  v  <  3.28W  ;  W  =  30  to  10  KT 

w 


(3.89) 


The  approximate  limits  of  X  on  the  pattern  center  line  (Y  =  0)  for  Equation 
3.88  is  XK  ±  r,  and  the  approximate  limits  of  X  in  Equation  3.89  is  X&  +  r  ±  r. 
The  hortest  arrival  time  is  for  particles  from  the  lowest  point  (center)  of 
the  cloud  bottom,  given  by  (h  -  b)/v5>  which  is 


t°  =  0.146W0'186  (1  -  0.0833W0136)  hr;  2  3.28W°'246 


(3.90) 


or 


t°  =  0.186W0,254  (4  -  0.0833W0-136)  hr;  v  <  3.28W0-246 
f  w 


(3.91) 


The  location  on  the  pattern  center  line  of  earliest  arrival  is  given  by 


V0  4.0 

X  =  v  t 
c  w  f 


(3.92) 


The  earliest  particle  fall  time  at  downwind  locations  outside  the  area  of 
the  limiting  ellipses  defined  by  Equations  3,98  and  3.89  is  "iven  by 


t,  -  X  /v 
f  e  w 


(3.93) 
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where  is  the  downwind  distance  along  the  pattern  center  line  to  the  center 
of  a  constant  arrival  time  ellipse.  The  distance  to  the  center  of  the  (half) 
ellipse  for  an  arbitrary  point  X,  Y  in  the  fallout  area  is  defined  by 

X  =  X  -  r%l  -  Y2/Y2  miles  (3.94) 

c  *  O 

In  general,  X  must  be  solved  by  an  iterative  process  using  Equations  3.94  and 
c 

3.85.  The  number  of  iterations  required  to  obtain  a  solution  of  X  and  Y  to 

co 

within  a  few  percent  of  the  converged  values  is  reduced  if  the  first  approxi¬ 
mation  of  in  Equation  3.85, is  taken  to  be  equal  to  X  -  r  if  Y  is  less  than 

r,  equal  to  X  if  Y  is  more  than  2r,  and  equal  to  X  -  0.5r  if  Y  is  hetween  r 
and  2r.  If  the  Iteration  is  done  by  computer,  the  last  given  first  approxima¬ 
tion  of  X  could  be  used  for  all  Y  values, 
c 

The  time  for  cessation  (or  last  time  of  arrival)  for  all  locations  in 
the  cloud  fallout  pattern  is  estimated  from  Equations  3.83  and  3.93  with  Xc 
estimated  from 

X  =  X  +  rWl  -  Y2/Y2  miles  (3.95) 

C  ”  o 

together  with  Equation  3.85.  If  the  center  point  of  the  constant  arrival  time 

ellipse  is  selected,  the  value  of  Y  from  Equation  '’.SS  is  obtained  directly 

o 

and  the  arrival  and  cessation  time  contours  are  then  calculated  from  Equations 
3.94  and  3.95. 

3.5.2  L  oendenc^  of  Particle  Fall  Times  and  Other  Parameters  on 
La  -eral  Wi  \d  Shear 

The  particie  fall  times  and  scaling  system  parameters  for  the 
larger  particles  tb<*t  fall  to  the  ground  from  stem  heights  are  assumed  to  be 
unaffected  by  lateral  wind  shear;  thus  lateral  wind  shear  effects  are  considered 
only  for  particles  falling  from  cloud  heights. 
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The  methods  described  below  for  estimating  the  particle  falling  speeds 
and  times  and  other  parameters  are  for  very  simple  lateral  wind  shear  patterns. 
The  basic  assumptions  for  developing  the  mathematical  representations  are: 

(1)  the  lateral  wind  shear  effect  occurs  only  through  the  height  range  of  the 
so-called  stabilized  cloud,  i.e.,  from  (h  -  b)  to  (h  +  b) ;  (2)  below  the 
height,  (h  -  b) ,  all  particles  fall  through  the  same  set  of  meteorological 
conditions  with  a  given  average  wind  speed  in  the  x  direction,  (thus  the  fall 
vectors  of  all  particle  groups  for  heights  less  than  the  bottom  of  the  cloud, 
are  parallel  to  the  x  axis);  (3)  the  lateral  wind  shear  through  the  height 
range  of  the  cloud  Is  either  described  by  a  given  angle,  6  ,  to  the  average 
direction  of  the  wind  below  the  cloud  height  or  by  a  single  valued  parameter, 
designated  as  <p,  in  degrees  per  unit  increase  in  height  above  the  bottom  of 
the  cloud;  (4)  the  magnitude  of  the  fall  vector  for  each  group  of  particles 
remains  constant  from  the  height  of  origin  in  the  cloud  to  the  ground;  and 
(5)  the  magnitude  of  the  wind  speed  is  constant  at  all  altitudes  (a  minor 
assumption  to  simplify  the  mathematical  representations). 

The  center  of  the  cloud  is  taken  as  the  center  of  the  coordinate  system; 
the  periphery  of  the  cloud  is  defined  then  by 

b(x+y)-r(b-z)=0  (3.96) 

The  time  of  fall  of  particles  with  a  given  value  of  v^  from  a  point  x, 
y,  z  j.n  the  cloud  to  tne  bottom  of  tr.e  cloud  is  given  by 

t.  -  (z  +  b)/v  (3.97) 

b  f 


or 

\  -  »>2  *  «» - 

v 

w 

where  X^,  is  the  location  of  the  particles  in  the  z  plane  at  (h  -  b)  at  the 
time,  t^.  The  lateral  displacements  are  thus  given  by 


,1/2 


(3.98) 
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X.  -  x  =  (v  /v  )(z  +  b)cos6 

D  W  f 


and 


Yb  -  y  =  (vw/vf)(z  +  b)sin8 


where 


6  =  (z  +  b)cp  degrees 


or.  If  8  is  constant, 


9  = 


9 

o 


degrees 


The  time  for  the  particle  group  to  fall  from  the  bottom  of  the 
the  ground  is  given  by 


t  *  (h  -  b)/v 
g  f 


or 


t 

g 


(X  -  X  )/v 
b  w 


where  X,  Y,  (Y  =  Y.  )  is  the  location  at  which  the  particles  arrived 
b 

ground  and  X  is  given  by 


X  =  X  +  (v  /v,)(h  -  b) 
b  w  f 

or,  in  combination  with  Equation  3.99,  by 

X  =  x  +  (v  /v  )[(h  -  b)  +  (z  +  b)cos91 
w  x 


(3.99) 

(3.100) 

(3.101) 

(3.102) 

cloud  to 

(3.103) 

(3.104) 

on  the 

(3.105) 

(3.106) 
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The  total  time  of  fall  is  equal  to 


+  t„- 
b  g 


given  by 


tf  =  (z  +  h)/vf 


(3.107) 


or 


t 


f 


r(x  -  x)2  t 

v 

w 


1/2 

<Y  -  y>] 


(3.108) 


It  should  be  noted  that  combinations  of  the  above  equations  do  not  provide 
explicit  solutions  for  cloud  periphery  values  of  x,  y,  and  of  z  in  terms  of 
and  the  cloud  geometry  equation  constants  when  S  depends  on  z  as  given  by 
Equation  3.101.  A  solution  for  v  ,  however,  is  possible.  (It  may  be  noted 
that  particles  with  the  same  value  of  originating  at  different  heights  in 
the  cloud  do  not  have  the  same  diameter;  for  particles  with  a  given  diameter, 
vf  is  also  a  function  of  z.) 

Estimates  of  the  fall  velocity  of  particles  originating  from  the  top  and 
bottom  of  the  cloud  volume  and  landing  at  the  point  X,  Y  may  be  made  from 


|  [  (h  -  b)X  +  (z  +  b)  (Xcos9  +  YsinQ )  ±^r2(l  -  zVb2)^  -  i|2] 


2  2 
X  +  Y 


2,,  2. .2. 

r  (1  -  z  /b  ) 


(3.109) 


where 


Vj  -  [(h  -  b)  +  (z  +  b)cos9]2  +  (z  +  b)2sin20 


(3.110) 


and 


\p2  =  (h  -  b)Y  -  (z  +  b)  (XsinS  -  YcosS)  (3.111) 

for  z  values  between  -b  and  b.  The  limiting,  or  boundary,  values  among  X,  Y, 
and  z  at  which  real  (as  opposed  to  Imaginary)  values  of  v  exist  are  obtained 
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by  setting  the  quantities  within  the  radical  of  Equation  3.109  equal  to  zero. 
The  deposition  pattern  boundary  lines  for  particle  groups  originating  at  the 
height  z  +  h  at  the  periphery  of  the  cloud  ellipsoid  are  thus  given  by 


Y 

m 


(z  +  b)Xsin9  ±  r  -  z2/b2)^ 
(h  -  b)  +(z  +  b)cos9 


(3.112) 


The  corresponding  values  of  the  fall  vector  may  be  estimated  from 


v  (m)  =  v 
f  w 


[ (h  -  b)  +  (z  +  b)cos9]  X  +  (z  +  b)Y  sin9 

m _ 

x2  +  Y2  -  r2(l  -  z2/b2) 
m 


(3.113) 


For  a  given  value  of  X,  Equation  3.112  represents  a  skewed  ellipt ically- 
shaped  curve  in  the  Y  ,  z  coordinate  system.  The  absolute  maximum  and  minimum 

R1 

values  of  Y^  depend  on  both  X  and  z.  Although  the  slope  of  the  lines  in  the 
X,  Y  plane  increases  as  z  increases  from  -b  to  b,  the  value  of  the  intercept 
at  X  equal  to  zero  has  a  maximum  for  particles  that  fall  from  the  center  height 
of  the  cloud  (z  =  O) .  Thus,  if  no  restriction  on  v^(m)  is  given  and  9  is  posi¬ 
tive,  the  largest  variation  would  be  2r  at  X  equal  to  zero  and  would 

m 

approach  the  value  2bXsin9/[(h  -  b)  +  2bcos9]  at  very  large  values  of  X. 

The  general  effect  of  lateral  wind  shear  on  the  boundaries  of  the  fallout 
patterns  and  on  arrival  and  cessation  times  can  be  investigated  through 
Equations  3.109  through  3.113.  Although  9  could  have  any  value  between  zero 
and  ±180  degrees,  the  average  lateral  wind  shear  at  heights  from  40,000  to  80,000 
feet  has  been  shown  to  be  correlated  with  wind  speed12  whereby  approximate  values 
of  Cp  can  be  estimated  from 

Cp  =  108  1,14  degrees/mile  (3.114) 

for  the  summer  season,  and 

-0  93 

cp  =  38.3  v  '  degrees/mile  (3.115) 

w 


46 


URS  702-1 


for  the  winter  season;  the  applicable  range  of  the  two  equations  is  for  z  +  h 
values  in  the  range  of  about  7  to  15  miles.  Thus,  for  an  integrated  wind  speed 
of  20  mi/hr,  the  average  lateral  shear  angle  from  the  bottom  to  the  top  of  the 
cloud  for  a  1  WT  yield  explosion  with  a  cloud  thickness  of  4.2  >>iles  would  be 
about  15°  for  the  summer  season  and  about  10°  for  the  winter  season.  The  data 
on  which  Equations  3.114  and  3.115  are  based  are  limited  and  perhaps  somewhat 
outdated;  further  studies  of  similar  data  would  be  desirable  to  determine  the 
frequency  of  occurrence  of  larger  lateral  wind  shear  angles  than  the  average 
angles  obtained  from  the  above  equations. 


If  9  is  assumed  or  taken  to  be  independent  of  z,  the  following  solutions 
of  z,  x,  and  y,  are  obtained  (the  solution  for  v^  is  the  same  as  that  given  by 
Equation  3.109  with  9  replaced  by  9^). 


vwb^  {  (Ysin9^  +  XcosO^'jv^  -  vw  Jb  +  (h  -  b)cos6^j| 


/  2  2 
(r  vf  ♦ 


2  2\ 
v  b  ) 
w  / 


± 


(3.116) 


r 
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where 


=  r2v2  +  v  b  |"(Ysin9  +  Xcos9  )v  -  v  (h  -  b)cos9  "1  (3.119) 

f  w  L  o  o  f  w  oJ 


and 


Y4 


.VfX  * 


v  (h 
w 


2,,2 

V 


(3.120) 


The  limiting  lateral  extent  of  the  fallout  deposition,  from  the  above 

equations,  is  defined  by  the  condition  that  the  Az,  or  Ax  and  Ay  segments 

within  the  cloud  be  equal  to  zero  (which  is  the  same  as  the  limit  for  obtaining 

real  values  of  x,  y,  and  z;  the  limiting  value  of  Y  for  a  given  value  of  v^ 

and  X  (with  stated  values  of  v  and  9  )  is  given  by 

wo 


v  bsin9  |r2v2  +  v  bcos9  fxv  -  v  (h  -  b)"l 
w  o  [  f  w  oL  f  w _ 

=  /  2  2  O  2T\ 

v,{r  v  +  v  b  cos  9  ] 
f  \  f  w  oJ 


’22  2  2  2a  j'  r 

r  v,  +  v  b  cos  9  -  Xv.  -  v  t 

-  bfl  -  cos9  VI 

1*1 

f  w  0  l  f  w  L 

V  o/J 

u 

(r2v 


2  2  2  2. 
.  +  v  b  cos  9 
f  w 


(3.121) 


An  explicit  equation  for  the  limiting  values  of  for  the  fallout  at  a 

given  location  is  not  possible  since  the  limiting  equation  contains  terms  of 

v^  with  exponents  from  zero  to  4.  The  values  of  v  that  give  maximum  and 

minimum  values  of  Y  at  any  value  of  X  may  be  estimated  from 
m 


v.(max)  =  (v  /X) 
f  w 


r  r2(h  -  b)b 


(h  -  b)  + 


+  bcos9  ) 
o 


2  I 


2r2(h  -  b)b2cos9  +  (X  +  bcos9  )2 

o  o  -J 


2bcos9  }  (3.122) 
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and 


vf(min)  =  (v^/X)  j(h  -  b)  +  — 


2r2(h  -  b)b2cos9 


r  (h  -  b)b  +  X  +  2(h  -  b)cos9 


J2) 


(3.123) 


Substitution  of  Equations  3.122  and  3.123  into  Equation  3.121  gives  the 
following  general  equation  for  the  boundaries  of  the  fallout  as  a  function  of 
the  downwind  (Us  ance,  X, 


bXsinS 

o 

1  2,2  | 
r  *.  *  1 

tv  - <h  -  »>] 

I  bXcos9 

o 

/  2,2 
\  vm 

+  b2X2cos29  ) 
o) 

1 

V('“  -  -2>2) 


r2.2  +  bVcos9 

-  x2| 

[*  -  r 

h  -  bfl 

.  -  cos9  Y|] 

l2l 

L  m 

o  1 

|  m  L 

V 

o/J  j 

u 

7  2,2  2  2  7  \ 

(r  1  +  o  X  cos9  ) 

\  m  <V 


(3.124) 


in  which 


=  Xv  (max)/v 
max  f  w 


(3.125) 


and 


m 


=  Xv  (min)/v 
min  f  w 


(3.126) 


as  obtained  from  Equations  3.122  and  3.123  respectively, 
and  (max)  refer  to  the  v  values  giving  the  smallest  and 
rather  than  the  smallest  and  largest  values  of  v^.. 


The  notations  (min) 
largest  values  of  Y 

m 


The  effect  of  lateral  shear  within  the  height  range  of  the  cloud  thickness 

is  illustrated  by  the  following  values  of  Y°  calculated  for  a  i  MT  vield  clouc 
o  m 

and  9^  value  of  30  (note  tnat  v^  does  not  appear  in  Equation  3.124): 
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X 

(mile3) 

y° 

max 

(miles) 

yo 

min 

(miles) 

AY° 

m 

(miles) 

AY°/2r 

m 

10 

10.3 

-8.1 

18.4 

1.01 

20 

11.5 

-7.2 

18.7 

1.03 

50 

15.5 

-4.9 

20.4 

1.12 

100 

23.2 

-3.0 

26.2 

1.44 

500 

94.0 

-0.1 

94.1 

5.17 

The  overall  pattern  spread  at  100  miles  away  from  ground  zero  for  30°  of  shear 
results  in  only  a  44  percent  Increase  in  the  apparent  pattern  width  relative  to 
the  "stabilized"  cloud  diameter.  Although  the  center  line  of  the  pattern  is 
slightly  curved,  the  angle  between  the  center  line  and  the  X  axis  to  a  distance 
of  50  or  100  miles  for  the  illustration,  is  about  6°. 

3.5.3  Dependence  of  Particle  Fall  Times  and  Their  Parameters  on 
Total  Wind  Shear 


For  overall  wind  shear  effects,  variation  in  wind  speed  as  well  as 

direction  over  the  range  of  cloud  thickness  must  be  taken  into  account.  To 

illustrate,  the  very  simple  shear  condition  is  assumed  in  which  v  is  the  inte- 

w 

grated  wind  speed  up  to  the  height  of  the  cloud  bottom  and  in  which  the  integrated 
wind  speed  increases  at  a  constant  rate  of  S  miles/hr  per  mile  at  greater  heights 
(a  variation  in  speed  with  height  may  also  occur  at  lower  heights,  which  would 
be  accounted  for  in  the  integrated  speed  to  the  cloud  bottom;  this  integrated 
wind  speed  would  be  approximately  the  same  for  all  v^  if  time  variations  in 
the  wind  speeds  and  directions  at  all  z  are  not  considered).  The  angle  between 
the  direction  of  the  integrated  wind  speed  vector  to  the  bottom  of  the  cloud 
and  the  direction  of  the  wind  above  the  bottom  of  the  cloud  is  designated  as  6. 

The  direction  of  the  integrated  wind  speed  vector  to  the  bottom  of  the  cloud  is 
defined  to  be  parallel  to  the  X  axis.  With  these  definitions  and  assumptions, 
the  wind  speed  over  the  height  range  of  the  cloud  is  represented  by 


v  (z)  =  v  +  (z  +  b)S 

w  w 


(3.127) 
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With  Equation  3.127,  an  explicit  solution  for  x,  y,  and  z  cannot  be 
obtained  for  the  given  cloud  geometry;  however,  the  values  of  v  for  particle 
groups  originating  at  the  height  z  +  h  on  the  cloud  periphery  and  landing  at 
the  point  X,  Y  may  be  estimated  from 


v  (h  -  b)X  +  fv  i-  (z  +  b)S~|  (z  +  b)(Xcos9  +  Ycos9) 

w  L  w _  J _ 


2  2  2 
X  +  Y  -  r  (1 


2  2 
z  /b  ) 


V- 


r2(l  -  z2/b2)t 


X2  +  Y2  -  r2(l 


2  2 
z  /b  ) 


(3.128) 


where 


=  |  v^(h  -  b)  +  £v^  +  (z  +  b)sj  (z  ■)  b)cos&J' 
+  £v^  +  (z  +  b)sj2(z  +  b)”sin20 


(3.129) 


and 


t  =  v  (h  -  b)Y  -  fv  +  (z  +  b)s"l(z  +  b)(Xsin0  -  Ycos9)  (3.130) 

6  w  L  w  J 

The  maximum  values  of  Y  for  ihe  particles  originating  at  the  height  (z  +  h)  and 
landing  at  the  downwind  distance  X  are  given  by 


Y 

m 


+  (z  +  b)sj(z  +  b)Xsin9  ±  r  -  z2/ >2)v, 
v^Ch  -  b)  +  [v  +  (z  +  b)S^(z  +  b)eos0 


(3.131) 


Although  an  explicit  equation  for  the  maximum  and  minimum  values  of  Y  is 

m 

not  readily  obtained,  approximate  values  can  be  derived  graphically.  The  effect 
of  total  shear  on  the  spread  of  the  pattern  is  illustrated  by  the  following 
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values  of  \'L  calculated  for  a  1  NfT  yield  cloud,  a  9  value  of  30°,  a  v  of 
m  J  o  w 

15  mi./hr,  and  an  S  value  of  2.1  per  hour: 


X 

(miles ) 

Y° 

max 

(miles ) 

Y 

min 

(miles ) 

oY° 

m 

(miles) 

<jY°/2r 

m 

10 

10.6 

-7  .9 

18.5 

1  .02 

20 

12.3 

-6.8 

19.1 

1.05 

50 

17.8 

-4  .6 

22.4 

1 .23 

100 

28.4 

-2.8 

31.2 

1.71 

500 

124.7 

-0.1 

124.8 

6.85 

The  overall  pattern  spread  at  100  miles  away  from  ground  zero  for  30°  of 
lateral  shear  and  an  S  value  of  2.1  ,ier  hour  of  height  above  the  bottom  of  the 
cloud  results  in  a  171  percent  increase  in  the  apparent  pattern  width  relative 
to  the  "stabilized"  cloud  diameter.  The  average  angle  between  the  center  line 
of  the  pattern  and  the  X  axis  at  a  downwind  distance  of  100  miles  is  slightly 
larger  than  7°.  The  values  of  i f  for  the  particle  group  that  land  at  X,  Y^ 
are  given  by  the  first  term  of  Equation  3.128. 

The  lateral  displacement  of  particles  with  a  given  value  of  v^.  originating 
at  the  height  z  +  h  on  the  cloud  periphery  and  landing  at  the  downwind  distance 
X  may  be  illustrated  from 


The  limiting  value  of  v^  for  ail  real  values  of  Y  in  terms  of  z  and  X  are 
given  by 
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v  (li  -  b)  +  fv  t  lz  t-  b)S*](z  -t  b)cos9 

vt(x>  ^  — - - J - 

*  «/  2  2 
X  ±  rVl  -  z  /b 


(3.133) 


And  the  corresponding  lateral  d  ispl  a<  cr.eni  limits  i  rom  Equation  3.132,  are 


fv  +  (z  +  b)sl  (  v  r  b)sin9  ,X  ±  r  W1  -  z“/bi| 

v  (x )  =  -L-— - i - k, - 1 _ ... — J 

.^(h  -  b)  +  f-.  i (z  +  b)S^(z  +  b)  eos9 


(3.134) 


The  equations  for  v  (X)  and  V(X)  give  the  envelopes  of  these  two  variables 

at  a  given  value  of  X  for  all  particle  groups  within  the  cloud  volume.  The 

paired  values  of  each  are  for  the  particles  originating  from  the  cloud  ellipse 

at  y  equal  to  zero.  The  general  shapes  of  the  v  envelopes  are  shown  in 

Figure  3.4  for  X  equal  to  50  and  100  miles.  In  addition,  specific  variations 

of  with  z  +  b  are  shown  for  the  X,  Y  point.-,  of  50,10  and  100,15.  For  each 

of  these  curves,  the  values  of  z  +  b  and  v,  that  give  the  minimum  and  maximum 

f 

of  v^.  and  t  are  indicated;  the  points  of  tangency  with  the  envelopes  are  also 
shown  as  well  as  the  major  parameter  values  used  ir.  the  calculations. 

An  explicit  general  solution  for  z  in  terms  of  other  parameters  is  not 
possible;  the  general  equation  of  z  is 

b2sV‘  +  2 bS ( v  b  +  2S)z3 

V 

■  i  *  O  '  ^  |  c 

.•  .  r“  -  2b_S(Xcoso  +  Ysin9)|+  v  b”  |v  +  2bv..  j)cosS  I  +  6b'JS(v  +  bS)}zil 

|t|_f  Jw[_w  J  wj 

+  2b2(v  +  2bS)|b(v  +  bS)  +  v  (h  -  b)cos9  -  v, (Xcos9  +  Ysin9)jz 

w  L  w  w  f  J 


+  b  ( v  +  bS )  b  ( v  + 
w  w 


|b(v^  +  bS)  -  2  jv ^.(Xcos9  +  Ysin9)  -  v^(h  -  b)cos9jj 
+  b2[xvf  -  vw(h  -  b)J2  +  v2fY2  -  r2)  =  0 


(3.135) 


Approximate  representations  of  the  fallout  pattern  boundaries  and  the 
general  directions  of  the  boundaries,  range  of  v  values,  and  times  of  arrival 
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(z  +  fc)  (miles) 

Fig.  3.4.  Variation  of  with  z  +  b  Showing  Envelopes  at  X  Equal  to  50  and  100  Miles  and 
Specific  Variations  for  X,  Y  Locations  50,  10  and  100,  15 
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and  cessation  at  locations  ma>  be  given  in  relation  to  the  particle  tra¬ 
jectories  from  the  center  of  the  cloud  originating  at  a  given  height  along 
with  the  dimensions  of  the  circular  disc  cnciosj  .g  the  area  covered  by  the 
particles.  The  latter  is  given  by 


(X  -  X  )2  +  (Y  -  Y  )2  =  r2(l  -  z2/b2) 
c  c 


(3.136) 


where 


v  (h  -  b)  +  fv  +  (z  +  b)s"|(z  +  b)cos^ 
w  L  w  _ J _ 


(3.137) 


and 

fv  +  (z  t  b)sl(z  +  b)sin9 

Y  =  ^ — - * -  (3.138/ 

c  vf 

The  paired  values  of  X  and  Y  ,  or  the  line  on  the  X,  Y  plane  along  which  the 

c  c 

particles  with  any  value  (originating  at  the  height  z  +  b  at  the  center  of  the 
cloud)  arc  deposited,  is  given  by 


Y 

c 


[vw  +  (z  +  b)sj(z  +  b)X^sin9 
v^(h  -  b)  -*■  [v^  +  ( z  +  b)S^J  (z  +  b)eos8 


(3.139) 


The  values  of  z  for  particles  with  a  given  range  in  vf  values  landing  at 

the  location  X  ,  Y  ,  from  rearrangement  of  Equation  3.139  may  be  estimated  from 
c  c 


z 


2S 


4v  S(h  -  b)Y 
_ w _ c 

^  sTn5  -  Y  cos6 


.) 


(v  +  2bS ) 
w 


(3.140) 


In  Equations  3.137  to  3.139,  the  angle  9  may  be  a  function  of  z;  however,  the 
separation  of  the  two  variables  in  Equation  3.140  indicates  the  assumption 
that  9  is  independent  of  z.  The  value  of  for  the  particle  disc  centered 
at  X^ ,  Y  ,  and  origination  at  the  height  z  +  h  in  the  cloud  is  represented  by 
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v  (h  -  b)sin9 
_  w _ o 

f  X  sln9  -  y  cos9 
coco 


(3.141) 


The  value  of  Y  ,  from  combination  of  the  above  equations  with  Equation 

3.138,  in  terms  of  X  and  Y  is  given  by 
c  e 


Y  =  Y  i  (r/X  )%'X2  +  Y2 
me  cl"'  c 


(3.142) 


The  time  of  fall  from  the  height  z  +  h  of  the  particle  disc  centered  at 
X  ,  Y  ,  may  be  estimated  from  Equation  3,107  or  from 


(X  sin9  -  Y  sin9  ^  /„  4v  (h  -  b)Y 

t.  =  \C  °  — r-  2S(h  -  b)  +Vv2  +  j-  - - — C-  *■  'S  (3.143) 

f  2v  S(h  -  b)sm3  f  w  /X  sins  -  Y  coso  \ 

w  o  ’  \  c  o  c  o/ 


It  may  be  noted  that  while  Equation  3.136  defines  the  area  covered  by  a 

circular  disc  in  the  X,  Y  plane  centered  at  the  location  X  ,  Y  ,  it  defines  a 

c  c 

skewed  elliptically  shaped  curve  in  the  Xc ,  Y^  plane  having  a  pseudo-center  at 

the  location  X,  Y;  such  a  curve  represents  the  loci  of  the  particle  disc 

centers  on  the  X ,  Y^  plane  for  those  discs  whose  radius  is  Just  equal  to  the 

distance  between  the  disc  center  and  the  location  X,  Y.  Thus  the  particles 

with  the  largest  and  smallest  values  of  v  and  t^  that  deposit  at  the  location 

X,  Y  are  those  on  the  edge  of  four  different  particle  discs  that  originate  at 

four  different  heights  and  whose  centers  are  at  four  different  X ,  Y  locations. 

The  values  of  X  for  particles  landing  at  X,  Y  originating  at  the  height  z  +  h 
c 

(and  located  at  the  edge  of  the  particle  disc)  is  given  by 

(v  (h  -  b)  +  fv  +  (z  +  b)S"l(z  +•  bjeos^j 

x  =  Lz - Ljz - j - L  x 

c  *5 

|v^(h  -  b)X  +  £v^  +  (z  +  b)sj(z  +  b)(Xcos9  +  YsinS)J 


v^(h  -  b)  +  +  (z  +  b)sj(z  +  b)cossJWr2(l  -  z2/b2);5 


(3.144) 
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The  corresponding  values  of  Y  may  be  calculated  from  Equation  3.139. 

The  variation  of  Y  with  X  for  the  cited  conditions  for  the  X,  Y  locations 
c  c 

50,10  and  100,15  is  shown  in  Figure  3.5.  The  X  ,  Y'  locations  of  the  particle 

c  c 

discs  that  are  first  and  last  to  arrive  (designated  t  and  t  ,  respectively), 

a  c 

that  originate  at  the  cloud  mid-height  (z  =  0),  and  that  originate  from  the 

maximum  and  minimum  heights  in  the  cloud  (7,  and  Z  ,  ,  respectively)  are 

max  min 

indicated  on  each  of  the  two  curves.  The  X  ,  Y  locations  of  the  discs  giving 

c  c 

the  maximum  and  minimum  values  of  vf  and  t  are  obtained  graphically  from  a 
plot  of  t  as  a  function  of  v .  The  respective  values  of  these  parameters  are 
as  follows: 


X  =  50,  Y  =  10 


X 

c 

(miles) 

Y 

c 

(mi’es) 

Vf 

(mi/hr) 

lf 

(hr) 

(z  +  b) 
(ml les ) 

50.39 

2.74 

2.552 

3.37 

0.836 

'55.99 

4.37 

2.406 

3.73 

1.200 

58.51 

6.96 

2.508 

3.83 

1.850 

42.53 

8.50 

4.191 

2.63 

3.262 

44.91 

10.32 

4.309 

2.70 

3.850 

49.49 

12.09 

4.081 

2.92 

4.154 

X  =  100,  Y  =  15 

100.56 

7.07 

1.319 

6.71 

1.080 

104.98 

8.39 

1  .288 

6.98 

1.225 

108.55 

12.25 

1.334 

7.13 

1.750 

92.04 

17.13 

1.868 

5.77 

3.003 

95.23 

20.06 

1 .926 

5.83 

3.468 

98.62 

21.36 

1 .890 

6.00 

3.584 

No  explicit  equations  have  been  derived  for  estimating  the  maximum  and 
minimum  values  of  v^.,  t^,  and  z  for  the  particles  deposited  at  an  arbitrarily 
selected  X,  Y  location. 
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3. 5. 4  Wind  Speed  and  Shear  Definitions  for  the  Fallout  Pattern 
Scaling  Systems 

4 

As  previously  discussed,  the  fallout  pattern  scaling  system 
coefficients  were,  to  a  large  degree,  evaluated  from  data  on  the  fallout 
distributions  from  test  shots  Jangle  S  and  Castle  Bravo.  Relative  wind  data 
for  these  two  detonations  are  applied  here  to  illustrate  some  definitions  of 
wind  speed  and  shear  that  are  applicable  to  the  scaling  system. 

In  general,  the  integrated  wind  speed  refers  to  the  totul  displacement 
of  an  imaginary  particle  of  undefined  size  falling  with  constant  speed  from 
a  given  height  (or  altitude)  to  another  height  or  to  the  ground.  Thus,  an 
integrated  wind  speed  may  be  directly  obtained  from  balloon  sounding  data. 

Where  the  latter  are  given  as  a  function  of  azimuth  and  speed  for  given 
altitude  increments,  then  v°(Z),  the  integrated  wind  speed  from  the  ground 
surface  to  the  height,  Z,  is  calculated  from 

i  =  n 

Ev  wZ  sin9 . 

i  1  l 

v°(Z)  =  -  (3.145) 

W  Zsin90 

n 

■where  is  the  thickness  the  ith  altitude  increment,  is  the  wind  speed 

of  the  itn  altitude  increment,  Z  if  the  height  or  altitude  of  the  nth  increment, 

and  9°  is  obtained  from 
n 

i  =  n 

£  Vin9i 

tan9°  =  - - -  (3.146) 

n  i  =  n 

£  vicos9i 
i  =  1 

The  average  (unidirectional)  wind  speed  is,  in  similar  units,  defined  by 
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-o 


w 


i  -  n 


■iMi 


i  =  1 


Z 


(3.1)7) 


The  i  lit  eg  rat  i<()  (or  overage)  particle  displacement  speeds  for  particles 
falling  from  large  heights  ore  usually  less  than  the  integrated  wind  speeds 
since  solid  particles  (and  especially  those  with  large  diameters)  fall  faster 
at  the  higher  altitudes  than  at  the  lower  altitudes.  The  integrated  particle 
displacement  speeds  arc  given  by 


v  (Z) 
w 


1  =  n 


v  It ^sinQ ^ 


sin9 
i  n 


(3.148) 


where  At  is  the  time  of  fall  through  the  ith  altitude  increment,  t  is  the 
i  n 

total  time  of  fall  through  the  n  altitude  increments,  Rnd  5  is  obtained  from 


tan9 

n 


i  -  n 


viAt1sin0i 


v  At cos9 


(3.149) 


both  At  and  t  ,  for  a  given  altitude  increment  and  total  distance  of  fall, 
in 

vary  with  particle  diameter. 

If  the  wind  speeds  and  directions  remain  constant  over  the  period  of  fall 

of  all  the  particles  of  interest  from  all  Z,  then  a  single  value  of  v  (Z)  and 

v°  for  each  Z  is  obtained.  This  condition  may  be  approximately  met  for  clouds 
w 

from  small  yield  detonations  and  for  t  values  of  less  than  2  to  3  hours.  Eut 
J  n 

for  clouds  from  large  vieid  detonations  where  t  is  greater  than  2  to  3  hours 

■  n 

even  for  fairly  large  particles,  it  is  much  less  likely  that  the  wind  speeds 
and  directions  will  remain  constant  in  time  at  even  one  location;  they  generally 
are  also  not  likely  to  be  constant  over  the  deposition  area  because  of  the 
resulting  larger  displacement  distances. 
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Tin-  results  from  displacement  calculations  (details  to  be  given  in  a 
following  report)  for  the  Jungle  S  Shot  and  Cuslle  Bravo  Shot  clouds  are 
summarized  in  Tables  3.0  and  3.7,  respectively. 

~n 

The  computed  values  oi  v  ('/)  ,  v  (Z)  for  particle  diameters  from  75  to 

w  w 

— o 

•100  microns  and  v  arc  all  about  the  same  for  the  Jangle  S  Shot  cloud;  also.  . 

w  | 

the'  9  values  differ  verv  little  from  G  for  the  integrated  wine,  speed, 
n  n 

However,  because  of  the  higher  altitudes,  variation  of  the  wind  sp  eds  and 
azimuths  with  time  (and  space),  and  large  changes  in  angle  with  altitude, 
the  values  of  v  (Z) ,  v  (Z) ,  and  v  for  the  Castle  Bravo  Shot  show  much  larger 

WWW 

differences;  the  same  is  true  for  9  and  9°.  The  maximum  difference  In  9  is 

n  n  n 

8.3  for  the  particles  from  the  Jangle  S  cloud  and  2w  <  for  the  1  0-micron 
particles  from  the  Castle  Bravo  cloud. 

The  accumulated  (or  integrated)  shear  factors  relative  to  the  height  of 
the  bottom  of  the  cloud  and  the  angular  disp’ acoment s  for  particles  of  different 
diameters  are  summarized  in  Tables  3.H  and  3.9  for  the  Jangle  S  and  Castle  Bravo 
Shot  clouds,  respectively.  The  total  shear  factor  for  the  Jangle  S  Shot  cloud 
is  quite  large,  mainly  duo  to  the  increase  in  wind  speed  with  height  above  the 
cloud  bottom;  in  addition,  the  values  of  both  S  and  9  are  about  the  same  lor  I 

all  the  particle  diameters  from  75  to  100  microns.  The  values  of  S  and  9  for  j 

the  particles  with  100-  and  200-micron  diameters  from  the  Castle  Bravo  cloud  I 

] 

include  time  and  space  variation;  the  variation  in  9  with  height  is  the  major  j 

source  of  shear  since  the  displacement  of  the  particles  starting  at  the  lop  of  i 

t  lie  cloud  is  almost  180  from  the  directions  of  the  hodograph  of  particles  , 

falling  from  the  bottom  of  the  cloud  to  the  ground. 

The  above  described  angular  and  total  shear  effects  are  included  in  the 
fallout  pattern  scaling  functions.  This  fact  should  be  kept  in  mind  in  applying 
the  functions;  in  no  case  are  the  scaling  functions  directly  applicable  to  a 
no-wind  shear  condition  (t.e.,  the  patterns  always  will  contain  residual  effects 
of  this  shear,  Irrespective  of  the  values  of  v  used  in  the  calculations).  If  ! 

a  single  wind  speed  and  direction  is  assumed  for  all  altitudes  (a  condition  I 
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Table  3.6 

SUMMARY  OF  INTEGRATED  WIND  AND  PARTICLE  DISPLACEMENT 


SPEEDS  AND 

Win.1  or  Particle 
Displacement 

DIRECTIONS 

0 

FOR  THE 

0.11 

JANGLE  S 

Z  -  h 

b 

SHOT  CLOUD 

,  (1G3  ft) 

1.11 

1.48 

2.11 

— c 

)  - 

V 

(Z)  or  v 

(Z) ,  (mi/hr) 

W  ’ 

•V 

Average  Wind  Speed 

25.1 

25.4 

27.8 

28.7 

30.0 

Integrated  Wind  Speed 

24  !  8 

25.1 

27.4 

28.3 

29.5 

Particle  Displacement 

Speed 

d  =  75  microns 

24  .6 

24.9 

27.2 

2  S .  0 

29.3 

d  =  100  microns 

24.6 

24.9 

27  .2 

28.0 

29.2 

d  =  200  microns 

24.6 

24.9 

27.2 

28.0 

29.2 

d  =  400  microns 

24.5 

24.8 

27.1 

27.9 

29.1 

r\  O  p 

*3  or  9 

(degrees ) 

n  n 

- 

Integrated  Wind 

184.8 

185.1 

188.5 

189.6 

191  .2 

Particle  Displacement 

d  =  75  microns 

184 . 7 

180.0 

188.3 

189.4 

191.1 

d  -  100  microns 

184.7 

185.0 

188.2 

189.4 

191.0 

d  =  200  microns 

184 .6 

185.0 

188.2 

139.1 

191.0 

d  =  400  microns 

184  .6 

184 .9 

188.2 

189.3 

190.9 

2 .96 


31.5 

30.9 


30.6 

30.6 

30.6 

30.5 


193.3 


193.1 
193.1 
193.0 
193  .0 
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SUMMARY  OF  (RELATIVE)  INTEGRATED  r.'IND  AND  PARTICLE  DISPLACEMENT 
SPEEDS  AND  DIRECTIONS  FOR  THE  CASTLE  BRAVO  SHOT  CLOUD 
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Table  3.8 

ACCUMULATED  SHEAR  FACTORS  AND  ANGLES  FOR  THE 
JANGLE  S  SHOT  CLOUD 


Particle  Diameter 

z  -  hb. 

(103  ft) 

(microns) 

0 

1.11 

1  .48 

2.11 

2.96 

v  (Z), 
w 

(mi/hr) 

75 

<24.6)a 

41.8 

42.5 

43.3 

43.7 

100 

(24.6) 

41.7 

42.6 

43.3 

43.8 

200 

(24.6) 

41.9 

42.7 

43.4 

43.8 

400 

(24.5) 

40.8 

42.5 

43.2 

43.7 

6  ,  (degrees) 

75 

(184.7)* 

15.0 

15.9 

16.6 

18.0 

100 

(184.7) 

14.9 

15.8 

16.6 

18.0 

200 

(184 .6) 

14.8 

15.8 

16.6 

18.0 

400 

(184.6) 

14.7 

15.7 

16.7 

18.0 

S  , 

(hr)"1 

75 

- 

81.4 

63.8 

46.6 

34.1 

100 

- 

81.5 

64.1 

46.7 

34.1 

200 

- 

82.3 

64.4 

46.9 

34  .2 

400 

• 

8°  .5 

64.6 

46.7 

34.2 

a  Values  in  parenthesis  are  for  the  integrated  particle  displacements 

f  rom  h.  . 
b 
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Table  3.9 

ACCUMULATED  SHEAR  FACTORS  AND  ANCLES  FOR  THE 
CASTLE  BRAVO  CLOUD 


z  -  hb 
;io3  ft) 

V 

(m 

(Z) 

w 

i/hr) 

6 

(degrees) 

S 

(hr'1 

) 

100  u 

200  u 

100  u 

200  u 

100  Li 

200  U 

0 

(19.6) 

3  ( 18 . 2 )a 

(254. 8)3 

(253. 0)a 

3.75 

23.7 

17.9 

70.2 

30.1 

5.73 

-0.46 

C  .75 

21 .9 

16.3 

64  .7 

69.6 

1 .35 

-1.16 

13.75 

22.4 

16.9 

72.6 

78.0 

1.07 

-0.51 

18.75 

17.9 

13.1 

86.8 

96.6 

-0.48 

-1.44 

23.75 

15.7 

12.8 

114.3 

131.0 

-0.88 

-1.22 

25.06 

16.4 

13.6 

116,9 

132.9 

-0.68 

-0.98 

28.75 

18.3 

15 . 6 

124.3 

138.2 

-0.25 

-0.44 

33.75 

19.4 

17.6 

138.1 

151.7 

-0.031 

-0,10 

38.75 

22.3 

20.9 

148.5 

159.4 

0.37 

0.36 

43.75 

24  .0 

23.3 

168.6 

178.0 

0.53 

0.62 

48.75 

27.9 

27.0 

169.5 

177.8 

0.89 

0.95 

50.12 

29.2 

28.0 

170.8 

179.0 

1.00 

1.03 

a  Values  in  parenthesis  are  for  the  integrated  particle  displacements 
from  fi 

D 
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that  corresponds  to  the  average  unidirectional  wind  speed  defined  above),  then 
the  computed  average  wind  speed  must  be  reduced  by  a  factor  of  15/28.7  (or 
0.52)  when  using  the  scaling  functions  for  yields  less  than  30  KT  and  by  a 

factor  of  15/26.6  (or  0.56)  when  using  the  scaling  functions  for  yields  larger 

* 

than  30  KT.  These  corrections  are  only  approximate  ones  for  the  larger  weapon 
yields;  the  correction  factor  actually  should  vary  with  particle  size.  However, 
these  corrections  will  give  a  first  order  approximation  of  the  pattern  scaling 
function  values  when  a  single  average  wind  speed  is  used.  Revisions  and  exten¬ 
sions  of  the  fallout  pattern  scaling  system  that  are  now  in  progress  should 
minimize  the  wind  speed  correction  problem;  the  major  question  to  be  resolved 
is  whether  the  adjustments  can  be  properly  introduced  without  also  introducing 
additional  and  more  involved  computational  steps  into  the  system. 

3.5.5  Distribution  of  the  Condensed  Radioactivity  Among  Particles  as 
a  Function  of  Falling  Velocity  Vector  of  Particle  Diameter 

The  activity  carried  by  particles  with  a  given  value  of  v  falling 
from  stem  altitudes  is  conveniently  represented  by 

i'O)  =  Kf(l)A^  R/hr  at  1  hr  (3.150) 

2 

where  A^.  is  the  activity  in  fissions/ft  for  a  unit  area  of  the  stem  (of  speci¬ 
fied  thickness)  at  the  altitude  from  which  the  particles  assume  gravity  fall, 

2 

and  Kf(l)  is  the  conversion  factor  from  fissions/ft  to  R/hr  at  1  hr  for  the 
particle  group.  The  total  activity  deposited  at  a  given  downvind  distance  in 
the  stem  fallout  pattern  is  calculated  from 

v (max) 

I  =  1.467v  I (K,(l)A;/v^)dv,  R/hr  at  1  hr  (3.151) 

s  w  1  f  f  f  f 

v^(min) 


*  No  correction  to  the  wind  speed  for  the  stem  fallout  pattern  is  suggested,  for 
most  applications  a  wind  speed  of  15  mi/hr  should  be  used. 
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where  v^(min)  and  v^.(max)  are  the  limiting  fall  velocity  vectors  in  ft/see 
for  the  smallest  and  largest  particles  that  are  deposited  (as  described  in 
Section  3,5.1),  and  v  is  the  average  v  nul  speed  in  mi/hr  from  the  surface 
to  Zj..  In  the  described  scaling  system  for  the  stem  fallout,  only  those 
particles  with  a  given  value  of  v  arrive  at  a  given  time;  the  arrival  time 
for  each  particle  group  Is  given  by  the  -sum  of  Equations  3.73  and  3.74  (cr 
3.75  and  3.76);  thus  the  activity  deposited  up  to  a  given  time  after  detonation 
is  calculated  by  graphical  integration  represented  by 


I 

s 


(t) 


v^(max) 


l.38v  I  (K  (l)A’.'v* 
w  I  f  f  f 


)Hv 


vf  (t ) 


R.  hr  at  1  hr 


(3.152) 


where  v^.(t)  is  the  fall  vector  of  the  particle  group  arriving  at  t  [v^(max) 
being  the  fall  vector  of  the  first  particles  that  arrive].  The  ionization 
rate  at  a  given  time  is  calculated  by  multiplying  Ig(t)  by  a  decay  correction 
factor,  df(t)  so  that 


I(t)  =  rif(t)Is<;)  (3.153) 

where  d^(t)  is  the  average  ratio  of  the  ionization  rate  at  the  time,  t,  to  the 
rate  at  H  +  1  for  all  the  particles  that  arc  deposited  at  the  location.  If  a 
large  variation  in  the  gross  decay  with  particle  sizes  is  to  be  assumed,  then 
d^.(t)  could  be  applied  in  Equation  3.152  prior  to  the  indicated  integrations. 

To  estimate  the  above  integrals  and  the  fallout  deposition  rates,  the 
variation  of  i^(l)  with,  v  must  first  be  computed  from  the  stem  fallout  pattern 
scaling  functions.  This  is  done  by  using  the  variation  of  I  with  X,  for 
X  2  (x^  +  Xg)/2  and  =  15  mph ,  along  the  center  line  of  the  stem  fallout 
pattern  and,  in  the  first  approximation,  taking  a  single  value  of  K^.(1)A^  in 
Equation  3.151  for  the  particles  that  fall  at  a  location.  The  single  (average) 
value  of  K^.(1)A^  is  then  written  as 
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I 

Kf(1)Af  =  - T - - - y  R/hr  at  1  hr  (3.154) 

1 ,467v  I l/v^(min)  -  l/v^(max) 

The  first  estimate  of  the  value  of  v^,  for  association  with  each  Kf(l)A^  value 
is  computed  from 

_  0  048  3  4 

log  vf  =  -0.050  +  0.769W  *  log  vf(min);  W  =“  10  to  5  x  lo4  KT  (3.155) 

or 


log  V  =  0.070  +  1.25w“0,043log  vf(max);  W  “  103  to  5  X  104  KT  (3.156) 

or  from  the  average  of  the  v  values  computed  from  both  equations.  The  adjusted 
2  1 

values  of  K^(l)A^/v^  as  a  function  of  v^  are  determined  from  successive  approx¬ 
imations  of  I  in  Equation  3.151  starting  with  K,(1)A'  for  K,(1)A'  and  v_  for  v  . 

s  f  f  f  f  f  f 

The  activity  deposited  by  particle  groups  with  a  given  fall  vector  from 
cloud  altitudes  at  a  given  downwind  location,  for  the  simple  fallout  pattern 
scaling  system  in  which  a  uniform  distribution  of  each  particle  throughout  the 
cloud  volume  is  assumed,  is  represented  by 


for  locations  on  the  fallout  center  line  with  either  a  directionally  balanced 
effective  (or  no)  lateral  wind  shear.  In  Equation  3.157,  A  is  the  average 

3  * 

initial  cloud  concentration  in  fisslons/ft  of  the  activity  carried  by  the 

particles  with  the  fall  vector  v  ,  K.(l)  is  the  conversion  factor  in  R/hr  at 
2  f  i 

1  hr  per  fission/ft  ,  and  the  remaining  combined  terms  represent  the  path 
length  through  the  elliptical  cloud  from  which  the  particles  with  the  fall 
vector  v^  originate  and  which  deposit  at  the  downwind  distance  X  on  the  pattern 
center  line.  Integration  of  1^(1)  over  all  the  particle  groups  that  land  at 
the  location  gives  the  total  activity  deposited;  this  is  represented  by 
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1 

s 


v  (max) 


1.467v 

w 


I 


(^.(D/v^ldv 

(min) 


R/hr  at  1  hr 


(3.158) 


for  v  in  mi/hr  and  v,  in  ft/scc. 
w  f 

The  average  deposition  rate  of  the  particle  group  with  the  fall  vector, 
v  ,  from  the  cloud  altitudes  is  given  by 

if(l)/Atf  =  3.6  X  103vfKf(l)Af  R/hr  at  1  hr/hr  (3.159) 

where  At^  is  the  time  period  over  which  particles  with  a  given  value  of  v  land 
on  the  ground;  £tf  may  be  obtained  from  a  plot  of  lf  as  a  function  of  v^;  for 
locations  near  the  pattern  center  line,  for  the  case  of  no  wind  shear,  it  is 
given  by 


At. 


'j[[ 


2btl(v2  +  v2b2r"2 
f  W 


)  -  ft  -  v.h*~’)x 

(v2  »  v2bV2) 

V  f  w  ) 


2  -2 
r 


(3.160) 


The  deposition  rate  of  all  the  particle  groups  arriving  at  a  given  time 
after  detonation  at  a  location  is  represented  by 


v  (max  ,  t) 


(dl  /dt)  =  5.28  X  10' 
s  t 


'«  f<K, 


(l)A^/v^.)dv^  R/hr  at  1  hr/hr 


(3.161) 


v^(nin , t) 


Integration  of  Equation  3.161  gives,  for  1^(0, 
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vf(max,t) 

(Kf(l)Af 

-  civ  R/hr  at  1  hr  (3.162) 

vf (mln,t) 

The  estimates  of  Kf(l)Af  for  difference  values  of  v  are  made  through  a 

series  of  approximations  starting  with  selected  values  of  1^  and  the  downwind 

distance  along  the  center  line  of  the  cloud  fallout  pattern  at  which  they  are 

calculated  to  occur,  using  the  simplified  scaling  system  with  a  15  mph  average 

wind  speed.  The  following  equations  may  be  utilized  to  provide  the  first 

estimates  of  Kf(l)Af  ar.d  vf  (designated  at  Kf(l)Af  and  vf,  respectively)  for 

each  set  of  paired  values  of  I  and  X. 

s 


I  (t)  =  5.28  x  10  v 
s  w 


L 

1 


Kf(l)Af 


4  0  411 

1.128  X  10  W  *  log  (p  (X) 


(3.163) 


where 


<Pf(X)  = 


(X  +  r)  W r2h2b~2  +  (X  +  r^ 
(X  -  r)  +Vr2h2b~2  +  (X  -  r)2 


,  X 


(3.164) 


and 


Cpf(X) 


(X  +  r)  -t 


X  + 


o 


Vr2h2b~2  +  (X  +  r)2 

Vr2h2b-2“7? 


X  S  r 


(3.165) 
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in  which 


r  -  0.464W0,431  miles,  W  =  1  to  10°  KT 

22-2  0  590  5 

r  h  b  =  31.00W  *  sq .  mi.,  W  =  30  to  10  KT 

X  =  0.186v  w0,254  miles,  W  =  30  to  105  KT 
o  w 

The  first  estimate  of  for  each  value  of  is  obtained  from 

—  70  0W°-164  5 

vf  =  — - — £ -  ft ''see  ;  X  >  r;  W  =  30  to  10  KT  (3.166) 

for  X  in  miles.  If  the  estimates  are  made  using  wind  speed  other  than  15  mi/hr, 

the  coefficient  70.0  is  replaced  with  4.67v  (for  v  in  mi/hr). 

w  w 

The  potential  exposure  dose  during  the  fallout  arrival  period  is  calculated 

by  integrating-  graphically  the  computed  values  of  I ( t )  from  t  to  t  for  the 

a  c 

fallout  from  both  stem  and  cloud  altitudes.  Also,  the  value  of  K^(l)  is 

adjusted  to  eliminate  the  instrument  response  factor,  l,  of  tquation  3.65; 

this  correction,  for  the  I  (t)  values  computed  through  use  of  the  simplified 

s 

fallout  scaling  system,  requires  increasing  th<  integrated  doses  by  a  factor 
of  1.33. 

Because  the  rather  complicated  form  of  the  above  equations  giving  the 
variation  of  Ig(t)  and  I(t)  with  time  after  detonation,  only  a  few  example 
calculations  were  made  from  which  approximating  functions  were  derived  to 
represent  the  variation  of  Ig(t)  or  the  accumulation  of  he  fallout  from  stem 
and  cloud  heights  with  time  after  detonation.  Since  the  model  of  the  stem 
depicts  the  fall  of  the  larger  particles  from  lower  heights  in  combination 
with  the  gradual  dilution  of  the  particle  concentrations  with  height  (and 
time)  and  the  deduced  activity-size  distribution  (from  the  arrangement  of  the 
stem  fallout  pattern),  a  fallout  rate  given  as  dlg(t)/dt  is  obtained  that  is 
highest  at  the  first  arrival  of  fallout  and  decreases  as  an  inverse  function 
of  time  during  the  fallout  period.  On  the  other  hand,  the  model  for  the 
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fallout  from  the  cloud  wherein  the  particles  are  assumed  to  fall  from  a 
stationary  source  volume  of  uniform  concentration,  the  fallout  rate  increases 
from  zero  at  the  first  arrlvul  of  fallout,  increases  to  a  maximum,  and  then 
decreases  to  zero  again  at  fallout  cessation.  This  general  picture  of  the 
deposition  rates  is  fairly  well  verified  by  analysis  of  observed  intensity, 
(exposure  rate)  data;  however,  detailed  data  from  time-interval  collectors 
show  large  short-time  variabilities  in  the  deposition  rates  which  indicates 
that  the  particles  are  not  uniformly  distributed  in  the  source  volumes. 

(Some  of  the  observed  variability  in  the  collector  data  is  undoubtedly  also 
caused  by  local  variations  in  the  surface  wind  speeds  and  directions  over  the 
deposition  period  and  the  apparent  efficiency  of  the  collectors.) 

The  simplified  versions  of  the  increase  of  I^(t)  with  time  during  the 

fallout  period  are  given  in  terms  of  the  fraction  of  1  thnt  is  deposited  at 

s 

a  given  fraction  of  the  deposit  period;  the  representation  is  defined  by 

I  (t)  =  f(t)I  (3.167) 

s  s 

and  the  fraction  of  the  time  during  the  deposition  period  is  defined  by 

T  =  (t  -  t  )/(t  -  t  )  (3.168) 

a  c  a 

For  the  iallout  from  stem  heights  with  f(t)  redesignated  as  f ^ ,  the  above- 
described  variation  of  the  fallout  rates  with  time  is  represented  approximately 
by 


f°Jl  -  (2t  -  l)2 

dfs/dr  =  -2-1 - - -  (3.169) 

o 

where  f  is  a  constant;  integration  of  Equation  3.169  under  the  condition  that 
f  is  zero  when  T  is  zero  and  f  is  equal  to  one  at  t  equal  to  one  gives 


f 

s 


(1/TT) 


-  (2t 


2  -1 
1)  +  sin  (2t 


1) 


0.500 


where 


f°  is  equal  to  1/tt. 


(3.170) 
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For  the  fallout  from  cloud  heights  with  f(t)  redesignated  as  f  ,  the 

c 

above-described  variation  of  the  fallout  rates  with  time  is  represented 
approximately  by 


df  /di 
c 


f ” V 1  -  <2-  -  l)2 


(3.171  ) 


for  a  symmetrical  rate  of  arrival  curve  that  has  a  peak  rate  of  urrival  at  the 
mid-time  of  fallout  period  (i.c.  ,  at  *  =  0.5).  Integration  oi  Equation  3.171 
gives 


=  (1/r)  (2t  -  1)^1  -  (2-  - 


l)2  t 


sin  1 ( 2 '  - 


1) 


0,500 


(3.172) 


In  regions  of  heavier  fallout  from  the  cloud  heights,  the  peak  fallout  rate 
generally  occurs  somewhat  earlier  in  the  fallout  period  than  at  the  mid-time, 
more  nearly  at  the  time  of  r  equal  to  0.4  than  0.5.  The  equations  for  f  with 
the  peak  rate  of  fallout  at  -  equal  to  0.4  are 


X 

c 


O/”) 


(2.5- 


V 


1)%1  -  (2. 


1) 


(2.5 


1) 


+  0.500;  -  S  0.4 


(3.173) 


and 

( 1 . 67 T  -  0 .67)  ^1  -  (1.67t  -  0 . 0 7 ) 2  +  sin_1(1.67-  -  0.67) 

+  0.500;  -  »  0.4  (3.174) 

The  calculated  variation  of  f  and  of  f  with  t  is  given  in  Figures  3.6 

s  c 

and  3.7,  respectively.  The  curves  can  be  used  directly  for  estimating  I(l); 

for  the  cloud  fallout,  and  the  values  of  f  taken  from  the  curves  for  T  equal 

c  o 

0.4  would  be  preferred.  At  the  shorter  downwind  distances  the  peak  fallout 
rate  may  occur  even  earlier  in  the  fallout  period;  also,  the  time  of  the  peak 


f  =  (I/-) 
c 
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ra^e  would  tend  to  shift  towards  the  mid-time  as  the  downwind  distance 
becomes  very  large.  The  approximations  for  f  and  f  given  by  Equations 
3.170,  3.173,  and  3.174  are  consistent  with  the  model  formulation  for 
application  to  computer  programming  in  damage  assessment  evaluations. 

To  illustrate  the  shapes  of  the  I(t)  curves  as  a  function  of  time  after 

detonation,  values  of  this  variable  were  calculated  for  the  following  assumed 

conditions;  for  the  stem  fallout,  t  is  0.5  hrs ,  t  is  1.0  hr,  I  is  1,000 

a  c  s 

R/hr  at  1  hr,  and  d(t)  is  equal  to  1/t ;  for  the  cloud  fallout,  t  is  1.0  hr, 

a 

t  is  2.0  hr,  I  is  1,000  R/hr  at  1  hr,  and  d(t)  is  equal  to  1/t.  Since,  in 
c  s 

general , 


t  =  t  +  (t  -  t  )T 

a  c  a 


C3.175) 


then 


t  -•  0.5  +  0.5T  (3.176) 

for  the  assumed  stem  fallout  case,  and 

t  =  1  +  T  (3.177) 

fer  the  assumed  cloud  fallout  case.  Also, 

I(t)  =  1,000  f  /t  (3.178) 


for  the  stem  fallout,  and 


I(t)  =  1,000  f  /t 
c 


(3.179) 


for  the  cloud  fallout.  The  calculated  values  of  I(t)  for  the  two  assumed  sets 
of  conditions  are  plotted  as  a  function  of  t  in  Figure  3.8.  The  peak  intensity 
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for  the  assumed  stem  fallout  case  Is  almost  1,100  R/hr  and  It  occurs  at  about 
0.8  hrs  (t  =  0.6).  The  peak  intensity  for  the  assumed  cloud  fallout  case  is 
505  R/hr  and  it  occurs  at  about  1.9  hrs  (r  =  0.9). 


The  air  ionization  rate  from  the  airborne  particles  can 

from  either  the  calculated  deposition  rates  or  the  estimates 

1 6 

From  calculations  made  by  Laumetz,  the  air  ionization  rate 
inside  a  plane  boundary  of  a  semi-infinite  volume  containing 
distributed  radiation  source  is  given  by 


be  estimated 
of  Kf(l)Af. 
at  three  feet 
a  uniformly- 


I  (t)  =  0.950  X  10_6A  A  (t)E  in  R/hr 
a  tot 


(3.180) 


in  which  At  is  the  activity  in  dis/sec  per  fission  at  the  time,  t,  after 

3 

fission,  A^(t)  is  the  concentration  in  fissions/ft  ,  and  E^  is  the  total 
gamma  ray  energy  in  Mev/dis. 

At  a  given  time,  At  and  E^  are  constant,  so  that  the  contributions  to 
activity  in  Af  from  all  the  particle  groups  arriving  at  a  given  time  can  be 
summed  as  in  the  case  of  the  deposited  material.  Multiplication  by  the  appro¬ 
priate  values  of  At  and  E^  for  the  designated  time  will  convert  the  sum  to 
R/hr.  The  fallout  from  stem  altitudes  arriving  at  a  given  instant  consists 
of  particles  with  a  single  value  of  v^  whose  air  activity  concentration  is 
estimated  from 


Ao(t)  =  KfdSvf(t)  (^d1)  fiSSl°ns/ft3  <3.181) 

The  values  of  A^  for  the  fallout  from  cloud  altitudes  arriving  at  a  given 
Instant  is  estimated  from 


A  (t) 
o 


3.38v 

w 


vf (max,t> 

^  (Af/Vf)(d  lo8  vf) 

vf (min.t) 


fissions/ft 


3 


(3.182) 
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in  which  Af  is  equal  to  if/Kf(l)  for  the  particles  arriving  at  the  time,  t, 
and  where  v  is  in  mi/hr  and  v  is  in  ft/sec. 

W  I 
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SUMMARY 


This  report  summarizes  the  fallout  pattern  scaling  relationships  that 
were  developed  in  the  period  1962  to  1964 ;  the  report  Includes  the  values 
of  the  scaling  equation  coefficients  that  were  derived  from  selected  fallout 
pattern  data.  The  meaning  of  the  scalar  wind  speed  multiplier  that  is  used 
in  the  scaling  equations  is  discussed  relative  to  computer  applications  of 
the  scaling  system  and  approximate  wind  speed  adjustment  factors  for  use 
with  wind  speed  averages  that  may  be  assumed  in  such  applications  are  pro¬ 
vided.  The  relative  degree  of  wind  shear  inherent  in  the  scaling  system 
parameters  is  also  discussed  in  some  detail.  Basic  equations  for  relating 
surface  density  of  radionuclides  and  air  ionization  rates  including  consider¬ 
ation  of  fractionation,  surface  roughness,  and  instrument  response  are  given 
and  discussed  together  with  the  influence  of  these  factors  and  others  on  the 
limiting  values  of  K  factors  that  represent  the  relative  amount  of  the 
radioactive  sources  contained  within  the  deduced  area  covered  by  the  fallout 
patterns.  Scaling  equations  and  data  are  also  presented  for  use  in  esti¬ 
mating,  for  any  location  in  the  fallout  region,  the  time  of  fallout  arrival, 
the  time  of  fallout  cessation,  the  variation  of  the  exposure  rate  (i.e.,  air 
ionization  rate  in  roentgens  per  hour)  with  time  during  fallout  arrival,  and 
the  total  potential  exposure  from  the  time  of  fallout  arrival  to  selected 
later  times. 
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Appendix  A 

SUMMARY  OF  SELECTED  SCALING  SYSTEM  PARAMETERS 


Provisional  values  of  several  of  the  fallout  scaling  system  model 
parameters  for  weapon  yields  In  the  range  of  1  to  30  KT,  not  provided  in 
the  main  text  of  the  report,  are  as  follows: 


X2  ="  0.0139W0,494  -  10.34W‘°-1C7J  miles;  W  =  1  to  30  KT 

0  494  T  -O  1  S7l 

X3  0.0139W  *  vw  +  10.34W  miles;  W  *  1  to  30  KT 

X„  **■  0.0973V  W°,55°  miles;  W  -  1  to  30  KT 

4  w 

3  0  414 

Z  =  3.43  X  10  W  feel;  W  =  1  to  30  KT 

2  1 3 

I  =  3.483  X  104v  “mw  W1*122;  W  =  1  to  30  KT 
4  t  J  w 

m  =  0.793;  W  =  1  to  30  KT 
w 

X,  =  -0.385W0'496  miles;  (v  =  15  mi/hr);  W  =  1  to  30  KT 
1  w 

0  404 

X  *  X„  -  0.139W  -log  I  ;  W  =  1  to  30  KT 

A  <*»  6  I  J 

0.467 

X.  =  0.0556V  W  miles;  W  =  1  to  30  KT 

5  w 

=  0.0894v  W0-481  miles;  W  =  1  to  30  KT 

6  w 

X,  =  0.0919v  w0,586  miles;  W  =  1  to  30  KT 

7  w 

0.596 

XQ  =  0.128v  W  miles;  W  =  1  to  30  KT 

8  w 

X„  =  1.956v  W0,319  miles;  W  =  1  to  30  KT 

9  w 

v,.r/h  =  8.34W-0,03  ml/hr;  W  =  1  to  30  KT 
-0.050 

v.r/h  =  5.19W  mi/hr;  W  =  1  to  30  KT 

D 

v?r/h  =  5.04W-0’155  mi/hr;  W  =  1  to  30  KT 


(A-l) 

(A-2) 

(A-3) 

(A-4) 

(A-5) 

(A-6) 

(A-7) 

(A-8) 

(A-  9) 

(A-10) 

(A-ll) 

(A-12) 

(A-13) 

(A-14) 

(A-15) 

(A-16) 
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o  2  2  -2  o 

The  values  of  KA ,  n^,  (v^r  b  )  ^ ,  and  for  the  yield  range  of  1  to  30  KT 
are  the  same  as  those  given  in  Table  3.1  for  application  in  Equations  3.28 
and  3.30. 

Calculations  of  fallout  arrival  times  for  average  wind  speeds  of  10,  20, 
and  30  mi/hr  at  selected  downwind  (X)  and  crosswind  (Y)  distances  for  weapon 
yields  of  1,  3,  and  10  MT  are  summarized  ir.  Tables  1A ,  2A ,  and  3A,  respectively. 
The  wind  speed  corrections  discussed  in  paragraph  3.5.4  do  not  apply  to  the 
arrival  time  calculations.  Similar  calculations  on  the  fallout  cessation  times 
are  summarized  in  Tables  4A,  5A ,  and  6A . 


Calculations  of  the  exposure  dose  during  the  fallout  period  and  from 
fallout  cessation  to  H  +  36  at  selected  X.Y  locations  are  summarized  in  Table 
7A  in  terms  of  Exposure  Dose-Standard  Intensity  Ratios  for  a  detonation  yield 
of  1  MT  and  an  average  wind  speed  of  20  mi/hr.  The  calculations  show  that 
D(36)/Ig  is  essentially  independent  of  Y  at  a  given  value  of  X,  The  values  of 
D(36)/Is  are  represented  very  closely  by 

-3 

D(36)/Ig  =  3.16  e"9-79  X  10  X;  X  *  85  km,  All  Y  (A-17) 


and 

-2 

-2  72  x  in  x 

D(36)/Ig  =  1.73  e  ;  X  *  85  km,  All  Y  (A-18) 

tor  W  =  1,000  KT,  v^  =  20  mi/hr,  and  X  in  km.  The  additive  exposure  doses 
for  times  greater  than  H  +  36  may  be  obtained  from  the  so-called  dose-rate 
multipliers  (DRM)  given  in  Volume  II  of  Reference  4.  The  variation  of  I(t)/I 

s 

with  time  after  detonation  for  downwind  distances  of  10,  20,  50,  150,  250,  and 
350  km  is  illustrated  in  Figures  1-A  through  6-A . 
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Table  1A 

SUMMARY  OF  TIME  OF  ARRIVAL  OF  FALLOUT  FROM  CLOUD  HEIGHTS 
AT  SELECTED  X,Y  LOCATIONS  FOR  W  =  1,000  KT 
(Time  in  Hours  After  Detonation) 

±  Y  (km) 


(km) 

0 

5 

10 

20 

40 

60 

80 

100 

120 

1.  V 

w 

-  10  mi/hr 

10 

0.94 

0.94 

0.94 

0.96 

1.12 

20 

0.95 

0.95 

0.96 

0.97 

1.07 

30 

1  .04 

1 .04 

1.06 

1.12 

1.36 

50 

2.28 

2.84 

2.29 

2.32 

2.44 

2  .65 

150 

8.50 

8.50 

8.50 

8.51 

8.55 

8.61 

8.71 

8.87 

9.14 

250 

14.71 

14.71 

14.71 

14.72 

14 ,75 

14 . 80 

14.88 

15.00 

15.18 

350 

20.93 

20.93 

20.93 

20.94 

20,96 

21 .01 

21.09 

21.19 

21 .36 

2 .  v  =  20  mi/hr 
w 


10 

0.49 

0.49 

0.51 

20 

0.50 

0.50 

0.52 

30 

0.56 

0.58 

0.68 

0.81 

50 

1.18 

1 . 19 

1.21 

1 .30 

150 

4.29 

4.29 

4.30 

4  .32 

4.40 

4.60 

250 

7.40 

7.40 

7.40 

7.41 

7.46 

7.57 

350 

10.51 

10.51 

10.51 

10.52 

10.56 

10.64 

10.81 

450 

13.61 

13  .61 

13.62 

13.63 

13.66 

13.74 

13.88 

550 

16.72 

16.72 

16.72 

16.73 

16.77 

16.84 

16.97 

650 

19,83 

19.83 

19.83 

19.84 

19.88 

19.94 

20.08 

750 

22.94 

22.94 

22.93 

22.95 

22.98 

23.05 

23.18 

3.  v 

w 

-  30  mi 

/hr 

10 

20 

0.49 

0.50 

0.52 

30 

0.51 

0.52 

0.55 

50 

0.82 

0.83 

0.87 

150 

2.89 

2.89 

2.90 

2.93 

3.09 

250 

4.96 

4.96 

4.97 

4.98 

5.06 

350 

7  .02 

7.03 

7.04 

7.05 

7.10 

7.24 

450 

9.10 

9.10 

9.”. 

9.12 

9.16 

9.26 

550 

11,18 

11.18 

11.1 1 

11.19 

11.23 

11.32 

650 

13.25 

13.25 

13.25 

13,26 

13  .30 

13.38 

750 

15.32 

15.32 

15.32 

15.33 

15.37 

15.45 

850 

17.39 

17.39 

17.39 

17.40 

17.44 

17.52 

950 

19.46 

19.46 

19.46 

19.47 

19.51 

19.59 

1050 

21.53 

21.54 

21.54 

21.54 

21.58 

21.66 

1150 

23.61 

23.61 

23.61 

23.62 

23.65 

23.73 
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Table  2A 

SUMMARY  OF  TIME  OF  ARRIVAL  OF  FALLOUT  FROM  CLOUD  HEIGHTS 
AT  SELECTED  X,Y  LOCATIONS  FOR  W  =  5,000  KT 

Y  (km) 


X  (km) 

0 

5 

10 

20 

40 

60 

80 

100 

120 

1  .  v 

w 

=  10  ml 

/hr 

10 

1.16 

1.17 

1.17 

1.19 

1.27 

20 

1.15 

1.15 

1.15 

1.17 

1.23 

30 

1.20 

1.20 

1.20 

1.21 

1.28 

50 

1.72 

1.73 

1.74 

1.78 

1.95 

2.23 

150 

7.94 

7.94 

7.94 

7.95 

7.99 

8.07 

8.18 

8.33 

8.54 

250 

14.15 

14.15 

14.15 

14.16 

14.19 

14  .24 

14.31 

14.41 

14.55 

350 

20.37 

20.37 

20.37 

20,38 

20.40 

20.44 

20.51 

20.59 

20.71 

2.  v 

w 

=  20  mi/hr 

10 

_ 

20 

1.29 

1.31 

- 

30 

1.15 

1.16 

1.17 

1.22 

50 

1.21 

1.21 

1.22 

1.26 

150 

4  .01 

4.01 

4.01 

4.04 

4.14 

4.36 

250 

7.11 

7.12 

7.12 

7.13 

7.19 

7.30 

7.48 

350 

10.22 

10.22 

10.23 

10.24 

10.28 

10.36 

10.48 

10.72 

450 

13.33 

13.33 

13.33 

13.34 

13.38 

13.44 

13.55 

13.72 

550 

16.44 

16.44 

16.44 

16.45 

16.48 

16.54 

16.63 

16.79 

650 

19.54 

19.55 

19.55 

19.55 

19.59 

19.64 

19.73 

19.87 

750 

22.65 

22.65 

22.65 

22.66 

22.69 

22.75 

22.83 

22 . 97 

23.27 

3.  v 

w 

=  30  mi/hr 

10 

0.58 

0.59 

0.65 

20 

0.55 

0.57 

0.58 

30 

0.56 

0.56 

0.58 

50 

0.65 

0.66 

0.71 

150 

2.70 

2.70 

2.71 

2.75 

2.95 

250 

4.77 

4.77 

4.77 

4.79 

4.88 

5.09 

350 

6.84 

6.84 

6.84 

6.86 

6.92 

7.04 

450 

8.91 

8.91 

8.92 

8.93 

8.97 

9.06 

9.26 

550 

10.98 

10.98 

10.99 

11.00 

11  .04 

11.11 

11 .25 

650 

13.06 

13.06 

13.06 

13.07 

13.10 

13.17 

13.29 

750 

15.13 

15.13 

15.13 

15.14 

15.17 

15.23 

15.35 

850 

17.20 

17.20 

17.20 

17.21 

17,24 

17.30 

17.40 

950 

19.27 

19.27 

19.27 

19.28 

19.31 

19.37 

19.47 

1050 

21.34 

21.34 

21.34 

21.35 

21.38 

21.44 

21.54 

21.79 

1150 

23.41 

23.41 

23.42 

23.42 

23.45 

23.51 

23.60 

23.84 

A  .4 


140 


8.86 

14.73 

20.86 
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Table  3A 

SUMMARY  OF  TIME  OF  ARRIVAL  OF  FALLOUT  FROM  CLOUD  HEIGHTS 
AT  SELECTED  X,Y  LOCATIONS  FOR  W  =  10,000  KT 

Y  (km) 


X  (km) 

0 

5 

10 

20 

40 

60 

80 

100 

120 

140 

1.  V 

w 

=  10  mi/hr 

10 

1.46 

1 .46 

1.48 

1  .48 

1.54 

1.65 

20 

1.44 

1 .44 

1 .44 

1.45 

1 .50 

1.59 

30 

1.45 

1.45 

1 .45 

1 .46 

1.51 

1.59 

50 

1.57 

1.57 

1.58 

1.59 

1.64 

1.75 

150 

6.93 

6.93 

6.94 

6.95 

7.01 

7.11 

7.24 

7.43 

7.66 

7.96 

250 

13.15 

13.15 

13.15 

13.16 

13.19 

13.25 

13.33 

13.43 

13.57 

13.74 

350 

19.36 

19.36 

19.37 

19.37 

19.40 

19.44 

19,50 

19.59 

19.69 

19.82 

2.  v  = 

w 

=  20  ml/ 

hr 

10 

_ 

20 

1 .64 

1.65 

1.69 

30 

1.48 

1.48 

1 .49 

1.54 

50 

1.44 

1.44 

1 .45 

1.48 

1.65 

150 

3.50 

3.50 

3.51 

3.55 

3.70 

3.97 

250 

6.61 

6.61 

6 .61 

6.63 

6.71 

6.84 

7.04 

7.39 

350 

9.72 

9.72 

9.72 

9.73 

9.78 

9.87 

10.00 

10.20 

10.54 

450 

12.82 

12.82 

12.83 

12.84 

12.88 

12.94 

13.05 

13.20 

13.42 

550 

15.93 

15.93 

15.93 

15.94 

15.98 

16.03 

16.12 

16.25 

16.42 

16.71 

650 

19.04 

19.04 

19.04 

19.05 

19.08 

19.13 

19.21 

19.32 

19.48 

19.71 

750 

22.15 

22.15 

22.15 

22.15 

22.18 

22.23 

22.31 

22.41 

22.55 

22.75 

3.  v 

w 

=  30  ml/lir 

10 

_ 

_ 

20 

- 

- 

- 

30 

- 

- 

- 

50 

1.51 

1.52 

1.54 

150 

2.36 

2.36 

2 .38 

2.44 

2.71 

250 

4.43 

4.43 

4.44 

4.47 

4.59 

4.83 

350 

6.50 

6.50 

6.51 

6.52 

6.60 

6.74 

7  .00 

450 

8.57 

8.57 

8.58 

8.59 

8.64 

8.75 

8.91 

550 

10.64 

10.64 

10.65 

10.66 

10.70 

10.78 

10.91 

11.14 

650 

12.71 

12.72 

12.72 

12.73 

12.77 

12.83 

12.94 

13.12 

750 

14.79 

14.79 

14.79 

14.80 

14.83 

14.89 

14  .99 

15.14 

15.43 

850 

16.86 

16.86 

16.86 

16.87 

16.90 

16.96 

17 .04 

17.17 

17.40 

950 

18.93 

18.93 

18.93 

18.94 

18.97 

19.02 

19.10 

19.22 

19.42 

1050 

21.00 

21.00 

21.00 

21.01 

21 .04 

21.09 

21.17 

21.28 

21.46 

1150 

23.07 

23.07 

23.08 

23.08 

23.11 

23.16 

23.23 

23.34 

23.51 
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Table  4A 

SUMMARY  OF  TIME  OF  CESSATION  OF  FALLOUT  FROM  CLOUD  HEIGHTS 
AT  SELECTED  X,Y  LOCATIONS  FOR  W  =  1,000  KT 

Y  (km) 


(km) 

0 

5 

10 

20 

40 

60 

CO 

100 

120 

1  .  v 

w 

=  10  mi/hr 

10 

1.62 

1.61 

1.60 

1.56 

1.28 

20 

2.24 

2.24 

2.23 

2.20 

2  .05 

30 

2.86 

2.86 

2.85 

2 .83 

2.72 

50 

4.10 

4.10 

4.10 

4.98 

4.01 

3.87 

150 

10.32 

10.32 

10.32 

10.31 

10.27 

10.21 

10.12 

9.98 

9.71 

250 

16.53 

16.53 

16.53 

16.52 

16.50 

16.44 

16.37 

16.25 

16.08 

350 

22.75 

22.75 

22.75 

22.74 

22.71 

22.66 

22.59 

22.48 

22.32 

2.  v  =20  mi/hr 
w 


10 

0.85 

0.84 

0.80 

20 

1.16 

1.15 

1.13 

30 

1.47 

1.47 

1.45 

1.37 

50 

2.09 

2.09 

2.08 

2.03 

150 

5.20 

5.20 

5.20 

5.18 

5.11 

4.93 

250 

8.31 

8.31 

8.31 

8.29 

8,25 

8.15 

350 

11.42 

11,42 

11.41 

11.40 

11.36 

11.29 

11.12 

450 

14.52 

14.52 

14.52 

14.51 

14  .48 

14.40 

14.26 

550 

17.63 

17.61 

17.63 

17.62 

17.58 

17.51 

17.38 

650 

20.74 

20.74 

20.74 

20.73 

20.69 

20.62 

20.49 

750 

23 . 85 

23.85 

23.84 

23.84 

23.80 

23  .73 

23.60 

3.  v 

w 

=  30  mi/hr 

10 

_ 

_ 

20 

0.80 

0.79 

0.74 

30 

1.01 

1  .00 

0.97 

50 

1  .42 

1.42 

1.40 

150 

3.50 

3.49 

3.49 

3.47 

3.32 

250 

5.57 

5.57 

5.56 

5.55 

5.48 

350 

7.64 

7.64 

7.64 

7.62 

7.57 

7.44 

450 

9.71 

9.71 

9.71 

9.70 

9.65 

9.56 

550 

11.78 

11.78 

11.78 

11  .77 

11.73 

11.64 

650 

13.83 

13.85 

13.85 

13.84 

13.80 

13.72 

750 

15.93 

15.93 

15.92 

15.92 

15.88 

15.80 

850 

18.00 

18.00 

18.00 

17.99 

17.95 

17.87 

950 

20.07 

20.07 

20.07 

20.06 

20.02 

19.95 

1050 

22.14 

to 

to 

►-* 

22.14 

22 . 13 

22.09 

22.02 

1150 

24.21 

24.21 

24.21 

24.20 

24.17 

24.09 
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Table  5A 


SUMMARY 

OF  TIME  OF  CESSATION  OF  FALLOUT  FROM  CLOUD  HEIGHTS 

AT  SELECTED 

X,Y  LOCATIONS 

FOR  W  = 

3,000 

KT 

Y 

(km) 

(km) 

0 

5 

10 

20 

40 

60 

80 

100 

120 

140 

1 .  v 

w 

=  10  ml/hr 

10 

2.16 

2.16 

2.15 

2,11 

1 .93 

20 

2  78 

2.78 

2.77 

2.74 

2.61 

30 

o  .40 

3.40 

3.40 

3.37 

3.27 

50 

4.65 

4.65 

4 .54 

4.62 

4.55 

4.41 

150 

10.86 

10.86 

10.86 

10.85 

10.82 

10.76 

10.67 

10.54 

10.36 

10.07 

250 

17.08 

'7.08 

1  7.07 

17.07 

17.34 

16.99 

16.93 

16.83 

16.71 

16.54 

350 

23.29 

23.29 

23.29 

23.28 

23.26 

23.22 

23.15 

23.07 

22.96 

22.81 

2.  v 

w 

=  20  mi/hr 

10 

1.11 

_ 

20 

1.43 

1 .42 

- 

30 

1.74 

1 .73 

1,72 

1.C4 

50 

2.36 

2.36 

2.35 

2.30 

)  50 

5.47 

5.47 

5.46 

5.45 

5.37 

5.20 

250 

8.58 

8.58 

8.57 

8.56 

8.51 

8.42 

8.25 

350 

11.68 

11  .68 

11.68 

11.67 

11,63 

11.56 

11 .44 

11.21 

450 

14.79 

14.79 

14.79 

14.78 

14.75 

14.68 

14.58 

14.41 

550 

17.90 

17.90 

17.90 

17.89 

17.86 

17.80 

17.71 

17.56 

650 

21  .01 

21.01 

21.00 

21.00 

20.97 

20.91 

20.82 

20.68 

750 

24.11 

24 .11 

24.11 

24.10 

24.07 

24.02 

23.93 

23.80 

23.50 

3.  v 

« 

=  30  mi/hr 

10 

0.77 

0.75 

0.69 

20 

0.98 

0.97 

0.92 

30 

1.19 

1 .  18 

1.14 

50 

1.60 

1 . 59 

1.57 

150 

3.67 

3.67 

3.66 

3,63 

3.49 

250 

5.74 

5.74 

5.74 

5.72 

5.65 

5.46 

350 

7.81 

7.81 

7.81 

7.80 

7.75 

7.64 

450 

9.89 

9.89 

9.88 

9.87 

9.83 

9.75 

9.56 

550 

11.96 

11.96 

11.96 

11  .95 

11.91 

11.84 

11.70 

650 

14.03 

14.03 

14  .03 

14  .02 

13.98 

13.92 

13.80 

750 

16.10 

16.10 

16.10 

16.09 

16.06 

16.00 

15.89 

850 

18.17 

18.17 

18.17 

18.16 

18.13 

18.07 

17.97 

950 

20.24 

20.24 

20.24 

20.23 

20.20 

20.15 

20.05 

L050 

22.32 

22.32 

22.31 

22  .31 

22.28 

22.22 

22.12 

21.87 

L150 

24.39 

24.39 

24.39 

24.38 

24.35 

24.29 

24.20 

23.97 
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Table  6A 

SUMMARY  OF  TIME  OF  CESSATION  OF  FALLOUT  FROM  CLOUD  HEIGHTS 
AT  SELECTED  X,Y  LOCATIONS  FOR  W  =  10,000  KT 

Y  (km) 


X  (km) 

0 

5 

10 

20 

40 

60 

80 

100 

120 

140 

1 .  V 

w 

=  10  mi/hr 

10 

3.15 

3.14 

3.14 

3.11 

2.98 

2.71 

20 

3.77 

3.77 

3.76 

3.73 

3.63 

3.42 

30 

4.39 

4.39 

4.38 

4.36 

4.27 

4.10 

50 

5.63 

5.63 

5.63 

5.61 

5.54 

5.41 

150 

11.85 

11.85 

11.84 

11.84 

11.80 

11.74 

11.65 

11.52 

11.36 

11.13 

250 

18.06 

18.06 

18.06 

18.05 

18.03 

17,98 

17.92 

17.83 

17.71 

17.57 

350 

24.28 

24.28 

24.28 

24.27 

24.25 

24.21 

24.15 

24.08 

23.98 

23.86 

2.  v 

w 

=  20  mi/hr 

10 

_ 

_ 

20 

1.92 

1.91 

1.89 

30 

2.23 

2.22 

2.21 

2.14 

50 

2.85 

2.85 

2.83 

2.79 

2.52 

150 

5.96 

5.96 

5.95 

5.93 

5.84 

5.67 

250 

9.07 

9.06 

9.06 

9.05 

8.99 

8.89 

8.73 

8.44 

350 

12.17 

12.17 

12.17 

12 .16 

12.12 

12.04 

11.93 

11.76 

11  .45 

450 

15.28 

15.28 

15.28 

15.27 

15.23 

15.17 

15.08 

14  .94 

14.74 

550 

18.39 

18.30 

IS. 39 

18.38 

18.35 

18.29 

18.21 

18.09 

17.93 

17,65 

650 

21.50 

21 .49 

21.49 

21.19 

21.46 

21.41 

21.33 

21 .23 

21.08 

20.85 

750 

24.60 

24.60 

24.60 

24.59 

24.57 

24.52 

24.45 

24.35 

24.21 

24.01 

3.  v 

w 

=  30  mi/hr 

10 

_ 

_ 

20 

- 

- 

- 

30 

- 

- 

- 

50 

1.92 

1.91 

1 .90 

150 

3.99 

3.99 

3.98 

3.95 

3.79 

250 

6.07 

6.06 

6.06 

6.04 

5.95 

5.77 

350 

8.14 

8.14 

8.13 

8.12 

8.06 

7.94 

7.71 

450 

10.21 

10.21 

10.21 

10.19 

10.15 

10.06 

9.91 

550 

12.28 

12.28 

12.28 

12.27 

12.23 

12.16 

12.04 

11.83 

650 

14.35 

14.35 

14.35 

14.34 

14.31 

14.24 

14.14 

13.98 

750 

16.42 

16.42 

16.42 

16.41 

16.38 

16.32 

16.23 

16.09 

15.81 

850 

18.50 

18.50 

18.49 

18.49 

18.46 

18.40 

18.32 

18.19 

17.97 

950 

20.57 

20.57 

20.57 

20.56 

20.53 

20.48 

20.40 

20.28 

20.09 

1050 

22.64 

22.64 

22.64 

22.63 

22.60 

22.55 

22.48 

22.37 

22.20 

1150 

24.71 

24.71 

24.71 

24.70 

24.68 

24  .63 

24.56 

24.45 

24.28 
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Table  7A 

SUMMARY  OF  EXPOSURE  DOSE -STANDARD  INTENSITY  RATIOS 

UP  TO  H  +  36  FOR  SELECTED  X,Y  LOCATIONS 

<W  =  1,000  KT,  v  =20  mi/hr )a 
w 


(km) 

Dl/Is 

V1* 

D(36)/I 

s 

X  (km) 

D,/I 

1  s 

V1. 

D(36)/I 

s 

Y  = 

0  km 

Y  = 

40  km 

10 

0.331 

2.54 

2.87 

150 

0.051 

1.10 

1.15 

20 

0.438 

2  .21 

2.65 

250 

0.033 

0.806 

0.839 

30 

0.446 

1.99 

2.44 

350 

0.023 

0.620 

0.643 

50 

0.241 

1.70 

1 .94 

450 

0.018 

0.483 

0.501 

150 

0.065 

1.08 

1.15 

550 

0.014 

0.377 

0.391 

250 

0.038 

0.801 

0.839 

650 

0.012 

0.289 

0.301 

350 

0.026 

0.617 

0.643 

750 

0.010 

0.215 

0.225 

450 

0.020 

0.482 

0.502 

550 

0.016 

0.375 

0.391 

Y  = 

60  km 

650 

0.013 

0.288 

0.301 

750 

0.011 

0.214 

0.225 

150 

0.024 

1.12 

1.14 

250 

0.024 

0.813 

0.837 

Y  = 

10  km 

350 

0.019 

0.623 

0.642 

450 

0.015 

0.486 

0.501 

10 

0.279 

2.61 

2.89 

550 

0.012 

0.379 

0.391 

20 

0.410 

2.24 

2.65 

650 

0.010 

0.291 

0.301 

30 

0.364 

2.01 

2.37 

750 

0.009 

0.216 

0.225 

50 

0.230 

1.70 

1.93 

150 

0.064 

1.09 

1.15 

Y  = 

80  km 

250 

0.038 

0.803 

0.841 

350 

0.026 

0.618 

0.644 

350 

0.009 

0.632 

0.641 

450 

0.020 

0.482 

0.502 

450 

0.008 

0.492 

0.500 

550 

0.016 

0.376 

0.392 

550 

0.007 

0.383 

0.390 

650 

0.013 

0.288 

0.301 

650 

0.006 

0.294 

0.300 

750 

0.011 

0.214 

0.225 

750 

0.005 

0.219 

0.224 

Y  = 

20  km 

30 

0.262 

2.06 

2.32 

50 

0.192 

1.72 

1  .91 

150 

0.062 

1.09 

1.15 

250 

0.036 

0.803 

0.839 

350 

0.025 

0.618 

0.643 

450 

0.019 

0.482 

0.501 

550 

0.016 

0.376 

0.392 

650 

0.013 

0.288 

0.301 

750 

0.011 

0.214 

0.225 

a  D^  Is  the  exposure  dose  to  fallout  cessation  Y 
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Fig.  A6.  Variation  of  I(t)/I  with  Time  After  Detonation 
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